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ABSTRACT 



Tl-E PURPOSE OF THIS THESIS IS TO INVESTIGATE A METHOD OF 
IMPROVING DECISION MAKING RELATIVE TO THE PROBLEMS CREATED 
BY OIL SPILLAGE. MANY COUNTRIES AROUND THE WORLD, INCLUDING 
INDONESIA, ARE PLAGUED BY INCREASING POLLUTION FROM THESE 
SPILLS. 

THIS THESIS USES A SIMULATION TO CONSIDER THE SPREAD AND 
DAMAGE CAUSED BY OIL SPILLS USING DATA FROM SAN FRANCISCO 
BAY. A PROJECTION OF SOCIAL COSTS FROM THESE SPILLS HAS BEEN 
MADE. 

FORMULATION OF A METHODOLOGY FOR DERIVING THE SOCIAL COST 
OF OIL SPILLS IS A PREREQUISITE IN REACHING OPTIMAL, RATIO- 
NAL DECISIONS IN MANAGING OIL POLLUTION. SUCH DECISIONS MAY 
INCLUDE THE ESTABLISHMENT OF A FINE STRUCTURE, D£ TERM I NAT I CN 
OF THE REQUIRED LEVEL OF CLEAN-UP AND IDENTIFICATION OF 
SOCIALLY SIGNIFICANT SPILLS. 
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I.INTROOUCTIGN 



f 



M4N HAS BEEN POLLUTING THE WATERS OF THE WORLD FOR YEARS. 
UNTIL THE ECOLOGY MOVEMENT AND THE PRESENT RESOURCE CONSERV- 
ATION PROGRAMt LITTLE CONCERTED EFFORT HAS BEEN MADE TO 
REDUCE POLLUTION. 

SINCE 1970, THE U.S. ENVIRONMENTAL PROTECTION AGENCY HAS 
PLAYED A MAJOR ROLE IN ATTEMPTS TO REDUCE THE FREQUENCY 
OF OIL AND HAZARDOUS SUBSTANCE SPILLS AND TO MINIMIZE ENVI- 
RONMENTAL DAMAGE CAUSED BY THOSE SPILLS THAT DO OCCUR. 

IN ADDITION, THE U.S. COAST GUARD HAS BEEN INCREASING ITS 
EFFORTS IN DETECTION AND CLEANING OF HAZARDOUS SUBSTANCES 
SPILLED INTO BODIES OF WATER. 

OVER 13, COO SPILLS OF OIL (REF-3) AND HAZARDOUS SUBSTANCES 
OCCUR ANNUALLY. SPILLED INTO RIVERS, STREAMS, COASTAL WATERS 
ESTUARIES AND LAKES, OIL SPREAD IN A MATTER OF MINUTES BY 
THE FORCE OF CURRENT INDUCED BY THE W IND ,S ALT N I T Y AND TIDES. 
SPILLS NCT ONLY REPRESENT WASTED RESOURCES BUT CREATE SOCIAL 
COSTS TO THE SOCIETY NEARBY DIRECTLY AND INDIRECTLY. 

CIL POLLUTION IS THE ALMOST INEVITABLE CONSEQUENCE OF THE 
DEPENDENCE OF GROWING POPULATION ON AN INCREASINGLY OIL- 
BASED TECHNOLOGY. 

BECAUSE OF THE LARGE QUANTITIES OFTEN INVOLVED IN SPILLS, 

THE EFFECTS ARE NOT ALWAYS COMPARABLE TC THOSE CAUSED BY THE 
CHRONIC POLLUTION OF INDUSTRIAL AND MUNICIPAL DISCHARGES. 
SOME OF THE EFFECTS ARE OBVIOUS, SUCH AS POLLUTED BEACHES, 
RIVERS DCTTED WITH OIL SLICKS, DEAD BIRDS AND FISH . 

BUT THE ECOLOGICAL EFFECTS FROM SPILLS ARE NOT CONFINED TO 
THE IMMEDIATE OR OBVIOUS SINCE OVER A LONG PERIOD OIL SPILLS 
COULD CHANGE THE COMPOSITION OF AQUATIC COMMUNITIES OR DAMA- 
GE THE ABILITY OF THE SPECIES TO SURVIVE. 

THIS STUDY WILL ADDRESS THE PROBLEM OF INCORPORATING A SO- 
CIAL CCST FIGURE ON THE CONDITIONS SURROUNDING AN OIL SPILL. 
IT EMPHASIZES THE PROBLEM OF GIL POLLUTION IN THE SAN FRAN- 
CISCC BAY AREA AND ITS RELATED SOCIAL COST, A PREDICTION OF 
SOCIAL CCST HAS SEEN MADE USING COMPUTER SIMULATION WITH 
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PROBABLE, 



BUT ARTIFIAL, 



INPUT DATA. 
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II. FREQUENCY OF OIL SPILL 



A. GIL FCLLUTIGN IN THE WATERS OF THE UNITED STATES. 



THIS SECTION WILL STATISTICALLY SHOW GIL POLLUTION IN 
U.S. WATERS: TYPE OF DISCHARGE, LCCATICN OF THE DISCHARGE 
AND THE SOURCES QP THE DISCHARGE. DATA PRESENTED IN TABLE 
ONE AND FIGURES GME THROUGH THREE IS BASED ON THE U.S. COAST 
GUARC pollution P-iCIDENT REPORTING SYSTEM (PIRSi. 

MACRO DATA OF OIL POLLUTION TRENDS IS GIVE^ 3Y TABLS-l WHICH 
CONTAINS STATISTICS FOR CALENDAR YEaRS 1971 TO 1974. 

OIL POLLUTION INCIDENTS CCCURING IN 1973 AND 1974 6Y AREA 
AND LGCATILN aRE SHOWN IN FIGURE-1 AND FIGURE-2. 

THE LEADING SOURCES EACH YEAR IN THE TOTAL VOLUME QISCHARGEO 
ARE THOSE WHICH HANDLE LARGS VOLUME OF PETROLEUM PRODUCTS, 
SUCH AS TANKERS, REFINERIES ETC. 

FlGUPE-2 SHOWS THE RELATIVE VOLUME DISCHARGED INCLUDING 
ONLY DISCHARGES OVER 100,000 GALLONS (OR 573,500 LITERS) 

IN 1973 and 1974. i 
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POLLUTION TRENDS IN ALL U. S. WATERS 
CALENDAR YEARS 1971 -1974 
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POLLUTION INCIDENTS BY AREA 



GREAT 





LAKES 

% 



VOLUME 

fig-1 



Source : U.S. Coast Guard P.I.R.S, 



(Ref-2) 



POLLUTION INCIDFNTS BY LOCATION 




NUMBER 




fig- 2 

Source ; U.S. Coast Guard P.I.R.S. (Ref-2) 



SOURCES OF DISCHARGES >100 ,000 GALLONS 




Calendar Year 1973 




Coast Guard P.I.R.S. (Ref~2) 



Source 



U.S . 



B. -SPILL PROeABILITY 



Tr£ WAJCRITY GF SPILLS ARE QUITE S.MALL . HOwcVER. IT IS 
THE VOLUME OF LARGE SPILLS THAI HEAVILY IMFLUENCES THE SIZE 
OF AVERAGE SPILLAGE. 

RELIANCE CN THE ANNUAL AVERAGE VOLUME SPILLED OVER A 
PROJECTED TIME PERIOD CAN 3E QUITE MISLEADING. 

IN REALITY, "HE ENVIPONMENT AL IMPACT OF OIL SPILLS OSPENCS 
BOTH CN THE FREQUENCY AND SIZE OF SPILLS. 

IT ALSC CEPcNCS ON THE RATE C-F SPILL DISCHARGE RELATIVE TO 
THE ABILITY TO CLEAN UP THE SPILL. 

BOTH OEV4NNEY AND PAULSON (REF-2) CONCLUDE THAT THE OCCUR- 
ENCE OF A POLLUTION INCIDENT IS ESSENTIALLY A RANDO'^ PROCESS 
AND CAN BE DESCRIBED BY A POISSON DISTRIBUTION : 



WHERE : 

N = # OF SPILLS 

X = VOLUME HANDLED 

\ = MEAN SPILL INCIDENCE RATE IN H S P IL LS / VOL UM E 

HANDLED. 

P(N|X) = PROBABILITY CF »N‘ SPILLS OCCUR IN3 GIVEN X. 



THEY CONCLUDE THAT : 

1. AVERAGE SPILL SIZES ARE RATHER MEANINGLESS STATISTICS 
SINCE THE VOLUME RANGE IS SO GREAT. 

2. THE TRUE IMPACT OF SPILLS IS A FUNCTION GF FREQUENCY 
SIZE AND LOCATION. 



BRUCE BEYAERT (REF-5) ESTIMATES THE RISK GF AN OIL SPILL BY 
MEANS OF A STATISTICAL PROBABILITY ANALYSIS USING THE PRO- 
CEDURE ILLUSTRATED Tm FTGURE-4, 

THIS ILLUSTRATION ASSUMES THAT AN ADEQUATE AND VALID BODY 
OF DATA IS AVATLiSLE TO INDICATE THE ACTUAL NUMBER AND SIZE 
OF ACCICENTAL GIL SPILLS . 

THE CATA WaS USED TQ DETERMINE BOTH THE PROBABILITY OF A 

THE PROBABLE DISTRIBUriGN GF SPILL SIZE 



SPILL EVENT AND ALSO 
FOR EACH EVENT. 

USING THIS ANALYSIS, THE 
SIZE CAN BE COMPUTED FOR 
THE RECURRENCE INTERVAL I 



RECURRENCE INTERVAL FOR ANY SPILL 
PANGS OF INTEREST. 

R the average period 



BETWEEN TWO SPILLS GREATER THAN OR EQUAL " 

THIS APPROACH IS EASILY UNDERSTOOD. BUT ITS _ _ _ 

THE VALID AND RELEVANT OIL SPILL DATA, WHICH UnFC RTUNAT ELY 
IS VERY HARD. TO DETAIN IN MOST INSTANCES, 



OF TIME 
4 SPECIFIED 
USE DEPENDS 



SIZE 

CN 
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ILLUSTRATION OF AN OIL SPILL RISK ANALYSIS 



1. SPILL frequency 

USING VALID AND RELEVANT HISTORICAL DATA, AN ESTIMATE FOR 
- N - IS MACE. THIS ESTIMA'^ICN IS THE AVERAGE NUMBER OF 
SPILL EVENTS PER YEAR: BASED ON CORRELATION WITH AN APPRC- 

. PRIATE OPERATING PARAMETER SUCH AS THE NUMBER OF TANKER PORT 
VISITS, THE NUMBER OF OIL TRANSFER OPERATIONS OR THE VOLUME 
OF OIL HANDLEC. 

2. SPILL SIZE DISTRIBUTION 

USING VALID AND RELEVANT HISTORICAL C ATA , OE T ERM INE THE 
PROBABILITY DISTRIBUTION OF ACTUAL SPILL SIZES. 




P* = PERCENT PROBABIL I TY, THAT A SPILL KILL HAVE A SIZE LESS 
THAN OR EQUAL TO S S GIVEN THAT A SPILL HAS CCCJRED. 

3. RECURRENCE INTERVAL 

FROM THE ABOVE, ■ R THE RECURRENCE INTERVAL IN YEARS CAN 
dE ESTIMATED FOR SPILLS OF VARIOUS SIZES BY USING THE FCLLCW 
ING RELATIONSHIP; 



R 



I 



100 

N ( 100 - s' ) 



FIGURE-4 



III. THE SPREAD OF OIL SPILLAGE IN WATER 



A. THE FATE OF OIL 



IF OIL IS SPILT ON LAND, PART OF IT WILL 3E A3SOR3cO INTO 
THE SOIL AND PART OF IT WILL FLOW OVER THE SURFACE SEEKING 
A LOW SPOT. 

THE FATE OF GIL SPILLED ON THE WATER IS A VERY COMPLEX 
SUBJECT. MUCH EFFORT HAS BEEN MADE TO ESTIMATE THE PHYSICAL 
SPREiDiNG AND MOVEMENT OF OIL ON THE SURFACE OF WATER UNDER 
THE INFLUENCE OF WIND, WAVES AND CURRENTS. 

THE PRECICTIVE MODELS DEVELOPED SO FAR ARE MOT CAPABLE OF 
HANDLING SUCH CCMPLEXlTIcS AS : 



1. EVAPORATION, 01 S SOLUTION, SEDIMENT AT ION, EMULSIFICATION, 
ALTO/PHOTO-OXIDATION AND 3 1 CDEGR ACAT I ON . 

2. THF change IN PHYSICAL AND CHEMICAL PROPERTIES OF 
FLOATING OIL RESULTING FROM PART ONE. 

3. THE EFFECT OF SEA CONDITIONS. - 



IT WAS OBVIOUS from EXPERIMENTS E EXPERIENCES THAT SMALL 
CUANTITIES OF OIL, I.E. A FEW TONS ( ONE TON OF OIL APPROX. 
EQUAL TC 6.5 BARRELS OR 250 GALLONS OR 947 LITERS) 

DISAPPEAR rapidly FROM THE MARINE ENVIRCNMENT. 

THE GENERAL NATURE OF GIL DISAPPEARANCE INVOLVES SUCH PhENC- 
MENA AS SPREADING, EVAPORATION, EMULSIFICATION, DISSOLUTION 
I ALTC-CXIlATIGN and 3I0DEGRA0ATICN (REF-10). 

IT GRADUALLY CISAFPcARS THROUGH DESTRUCTIVE AND DISPERSIVE 
PROCESS E5, LEAVING an ASPHALTIC MASS. 

OIL SINKS AFTLR ITS DENSITY IS INCREASED BY EVAPORATION, 

BY SCLUTICN OF ITS VOLATILE FRACTIONS, BY INCLUSION OF 
PARTICULATE MATERIAL AND BY OXIDATION (203ELL 1964, PIL?EL 
1968, REF-5). SINKING MAY HAVE 3EEN OF PARTICULAR IMPORTANCE 
IN THE SAN FRANCISCO SPILL (REF-5). 

A GREAT MANY PUBLICATIONS ON THE MStABCLISM OF HYCROCARBQNS 
IN WATER SUGGEST THAT ALL MOLECULES PRESENT IN CRUDE OIL 
CAN BE ATTACKED BY ENZYMES CF THE MI C RCGRGANI S MS . THIS 
NATURAL PROCESS , HOWEVER , MAY BE TOC SLOW UNDER NORMAL 
CONDITIONS AND THEREFORE , UNABLE TO CAUSE DESTRUCTION OF 
LARGE CIL SPILLS BEFORE DAMAGE IS DONE TO MARINE LIFE AND 
COASTAL AREAS. 

IN THE ALTO-aXIDATION PROCESS, TEMPERATURE IS AN IMPORTANT 
PHYSICAL FAC7GR. AT SEA, BELOW FIVE DEGREES CELCiUS, OXIDA- 
TION IS VERY SLOW, IT THUS OCCURS WITH GREAT DIFFICULTY IN 
LATITUDES ABOVE 75 DEGREES NORTH OR BELOW 75 DEGREES SOUTH. 
IN ECUATGRIAL REGIONS, THE RATE OF OXIDATION MAY ATTAIN 
A RATE CF SEVERAL HUNDRED GRAMS OF GIL PER CUBIC VETER OF 
HATER PER YEAR { G.D. FLOODGATE MI CRG E I AL DEGRADA . luN' CF 
OIL • .MAR. POLL BUL. 3 1972, P41-43 ) 

OIL REMAINING AT SEA FOR THREE MONTHS OR '^CRE, LOSES VOLUME 
continually and may be reduced TO asphaltic residle .REPRE- 
SENTING AS LITTLE AS 15? OF THE ORIGINAL VOLUME (SMITH 1968, 
REF-5 ). 

SIVACIEF AND MIKCLAj (REF-10) AFTER ASSESSING THE ROLE THAT 
VARIOUS FACTORS MIGHT PLAY IN DETERMINING THE FATE 0= OIL 
SPI LLS , ICENTI FI ED EVAPCRATICN AS BEING THE MOST SIGNIFICANT. 
EVAPORATION IS ENCHANCED BY INCREASING WIND SPEED, SEA SUR- 
FACE ROUGHNESS, AIR TEMPERATURE AND DECREASING Oil! FILM 
THICKNESS, lighter WEIGHT OIL EVAPORATES FASTER THAN HEAVIER 
TYPES. 



I 
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8. mcdel of disappearance of oil spills due to eva^craticn. 



TV<G 6ASIC ASSUMPTIONS MADE BY SIVADI5R AND MIKGLAJ IN 
THEIR MCDEL : 



1 iCiRT FPCiM SPREADING* THE ONLY PROCESS OCCURING TO 

-0 BIO- 

degradation are CONSICEKcD NEGLIGIBLE. 



2. OIL CONSISTS OF TWO ARBITRARY DEFINED PARTS i 

A. A VOLATILE FRACTION FROM WHICH ALL EVAPORATIVr 
LOSSES OCCUR 

B. A RESIDUUM FRACTION WHICH IS TOTALLY UNAFrECTED BY 
WEATHERING. 



F = Cl’S'T/d + C2-T) 



WHERE: 



c - TWP ■WFTGHT FRACTIF>N OF THE WEATHERED OIL 
■ SAMPlI WHICH IS EVAOGRATEO, IK =ERC:NT. 

n r? - CONSTAN'^t IN WHICH THE VALUE DEPENDS 
On'tHE TYPE OF OIL, WEATHER aMD WA.ER 
CONDITION. 

T = TIMEt IN MINUTES. 

AS TIME APPROACHES INFINITY, THE VALUE OF 'F' CONVERGES iQ 
20 - 22? (REF-10). 

SIN'*£ THE REMAINING OIL SLICK RESIDUE WOULD n^R iPL^ 

IrAVIT^ NEARLyShE SA.MB AS SEA !{|TER. THtRt I| C...N.IO 

PORT MECHANISM. 

C. MCDEL TO DETERMINE THE LEEWAY OF OIL oLiCKS. 

t^WAY^'^IS DfIiN£c''‘as‘tHE MOVEME^^ OIL SLICK OVER THE 
SaTER CLE TO THE ACTION OF THE WIND. 

SMITH (REF-4) CONCLUDE THAT : 

‘TMiliilalSiSliSBiSl 

ICN : 



OSL = Q.0179’!‘W10 + 0.0196 
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WhEPE 



OSL = GIL SLICK LEEWAY, IN KNOTS. 

WIO = WIND SPEED AT 10 METERS ELEVATION, IN KNOTS. 



3. THE EXPRESSION FOR WIND IN RANGE LESS THAN FIVE KNOTS 
ShCULD LSE: 



CSL = C.0199=i'W10. 



4. OIL SLICKS MOVE IN THE DIRECTION OF WIND ACROSS THE 
WATER surface. 

5. OIL SPILL VOLUME WAS NOT FOUND TO AFFECT THE MAGNITUDE 
OF THE SLICK LEEWAY, BUT VERY THIN OIL FILMS WERE 
FOUND TC EXHIBIT LITTLE OR NO LEEWAY. 

6. OIL SLICK LEEWAY INCREASES AS A POSITIVE FUNCTION 

OF SEA STATE, BUT THE RELATIONSHIP WAS NOT QLANTITATIVE 
LY DEFINED. 

7. GIL SLICK MOVES ACROSS THE WATER SURFACE UNCER THE 
INFLUENCE OF THE WIND LEAVING A THIN FILM ALONG THEIR 
PATH. 



C. OIL SLICK SPREADING ANALYSIS. 



IT IS A COMMON OBSERVATION THAT OIL, WHEN SPILLED ON WATER 
TENDS TC SPREAD OUTWARD ON THE WATER SURFACE IN THE FORM 
OF A THIN LAYER. THIS TENDENCY TO SPREAD IS the RESUL" O'F 
TWO PHYSICAL FORCES: THE FQPCE jF GRAVITY WHICH CAUSES THE 
LIGHTER CIL TC SEEK A CONSTANT LEVEL BY SPRE-AOING H- = IZONTAL 
LY,JLST AS IT WOULD ON A PLANE HORIZONTAL SOLiO SURFACE INC 
THE SURFACE TENSION FORCE OF PURE WATER, WHICH IS USUALLY 
GREATER THAN THAT OF THE OIL FILM FLOATING ON WAT?F. 

IN THE OPEN sea, THIS SPREADING TENDENCY IS AIDED BY WATER 
SURFACE MOTIONS INDUCED BY WAVES, WIND AND TIDAL CURRENTS. 
THESE EXTERNAL FORCES MAY BE THE ONLY CAUSE OF SPREAD IF 
THE SLICK HAS BEEN BROKEN INTO SMALL INCEPENDEN"^ PIECES. 

THE SPREAD RESULTING FROM THIS RANDOM MOTION Or THE SEA 
SURFACE IS VERY DIFFICULT TO ESTIMATE. MOST OILS SPILLED 
ARE MIXTURES CF COMPONENTS HAVING VARYING VAPOR PRESSURE 
AND SOLUBILITY IN WATER. THE LIGHTER MOLECULAR WEIGHT ARE 
MORE VOLATILE AND SGLUBLE,ANO THE HEAVI ER/RES I OUh IS DENSER 
AND MORE VISCOUS. 

FAY (REF-6,10) HAS DEVELOPED AN ANALYSIS FOR THE SPREADING 
OF A ONE DIMENSIONAL AND AXISYMETRIC OIL SLICK AS A FUNCTION 
OF TIME. THE ANALYSIS IS RESTRICTED TO A FIXED AMOUNT OF 
CIL IN THE initial SPILL AND TH^ WATER IS FREE OF MOTION 
INDUCEC EY WIND, WAVES AND TIDAL CURRENTS. 

THE SPREADING PROCESS PROCEEDS THROUGH THREE STAGES IN WHICH 
VARIOUS FORCE PAIRS ARE IMPORTANT AND THE SLICK FINALLY 
REACHES A TERMINAL SIZE. FOR THE AXI-SYMETRIC CASE, THE =A- 
DIUS CF THE SLICK 'R« IS RELATED TO TIME AFTER THE SPILL 
BEGINS , 'T* , BY THE FOLLOWING EQUATIONS: 
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GRAVITY-INERTIA STAGE; 



Rj = (A-GVT^) C, 



GRAVITY-VISCOUS STAGE: 



Rs= ( 



AG 

— ) . Q 



V 



to 



SURFACE TENSICN-ViSCGUS STAGE: 







1/4 




WHERE : 

Cl 

C2 

C2 

G 

A 

ao 

Po 

V 

a 

'"to 

T 

R 



1.14 

1.45 

2.30 

GRAVITY, GR/CM**2 
^i)” ^ Ri) 

WATER DENSITY, GR/CM^^J'S 
OIL DENSITY, GR/CN*»3 
VOLUMEt CM=i'*3 

SPREADING COEFFICIENT, OYNS/CM 
WATER KI NE NAT I C VISCOSITY, CM^*2/S5C. 
TIME, SECOND. 

RADIUS, CM 



THE FINAL RADIUS CF THE SPILL : 



R 



r 



n 



1/2 



WHERE: 

RF = RACIUS, METER. 
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V = VOLUME, 



THE TIME AT WHICH THE TRANSITION FROM THE INERTIA STAGE TO 
THE VISCOUS STAGE OCCURS ,• T 12 ' ,CAN BE FOUND BY EQUATING FROM 
THE SPILL RADII FROM EQUATION ’Rl» AND 'R2'. 





. / ^2 ■ 


) { ^ 




■ ■ 


' ^ A G 


V ^ 

(Jj 


SIMILARLY 


FOR THE 


VISCOUS- 


-SURFACE TENSICN TRANSITION TIME: 




IN THIS STUDY THE FOLLOWING VALUES ARE USED : 
G = 98C CM/SEC- . 

O' = 10 OYNE/CM 
P =1.0 

(jj 

pQ = 0.95 

LATITUDE = 38 DEGREE NORTH 
= 0.01 CM^ ,/SEC. 



C. MOVEMENT OF SPILLED GIL IN SAN FRANCISCO BAY AREA. 



THE INTENT OF THIS STUDY IS TO PRESENT DESCRIPTIONS OF 
GENERAL SURFACE WATER MOVEMENTS IN THE SAN FRANCISCO BAY 
AREA ANC ITS IMPACT ON THE DISTRIBUTION AND DISPERSAL OF OIL 
SPILLS IN THE BAY. 

A SIMULATION METHOD IS USED TO PREDICT THE DISPERSAL OF 
SPILLED OIL. 

THE SAN FRANCISCO BAY SYSTEM CAN BE DIVIDED INTO A N^RTH BAY 
- CENTRAL BAY AND A SOUTH BAY. 

PERRENIAL ESTUARINE CIRCULATION IN THE NORTH ANC CENTRAL BAY 
ANC THE ADJACENT OCEAN IS MAINTAINED BY THE SACRAMENTO - SAN 
JOAQUIN RIVER RUN OFF. ANNUALLY THE AVERAGE SURFACE SEAWARD 
DRIFT EXCEEDS FIVE KM PER DAY (REF-5). 

WIND IS A MAJCR FACTOR IN DETERMINING CIRCULATION IN THE 
SOUTH BAY. LOCAL WINDS ARE TOPOGRAPHICALLY CONTROLLED BY 
PREOCMINANTLY NORTHWEST WINDS DURING SUMMER AND THE SOUTH 
EAST WINDS DURING WINTER STORMS, WITH THE PREVAILING SUMMER 
WIND GENERALLY TWICE AS STRONG AS THAT IN WINTER (REF-5). 

THE SOUTH BAY CIRCULATION IS RELATIVELY SlUGGISH.THE ANNUAL 
AVERAGE SPEED OF SURFACE OR I FT , REGAP OL E SS OF DIRECTION IS 
BETWEEN CNE iiND TWO KM PER DAY. THE MOVEMENTS OF SURFACE 
DRIFTERS IN SAN FRANCISCO BAY ARE SHOWN IN FIGURc-5 (REF-5). 
AN EXAMPLE OF THE SPREAD AND THE MOVEMENT OF OIL SLICKS 
PLOTTED BY THE COMPUTER SIMULATION IS SHOWN IN FIGURE-6. 
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RNRLV3IS OF OIL SPILL MQUEMENT 
IN SRN FRRNCISCO BRV RRER 
fiS PREDICTED BV ESTURRINE NON TIDRL DRIFT 




IsEST SOUTHERN 




IV. FACTORS AFFECTED BY THE OIL POLLUTION 



CRUDE OIL AND OIL PRODUCTS SPILLED IN NATURE ARE ALTERED 
BY cVAPCRATION, BY DISSOLUTION, BY BACTERIAL AND CHEMICAL 
ATTACK. 

IN SPITE OF COMPLEX PROCESSES OCCURING CURING WEATHERING, 
MANY COMPOSITIONAL PARAMETERS ARE RELATIVELY STABLE AND ARE 
NOT OBLITERATED UNTIL AN ADVANCED STAGE OF DEGRADATION HAS 
BEEN REACHED (REF-1). 

THE STABLE PARAMETERS MAY AID THE IDENTIFICATION OF AN OIL 
POLLLTANT AND IN THE CORRELATION WITH ITS SOURCE FOR MANY 
WEEKS AFTER THE SPILL (REF-1, 2). 

THE EFFECTS OF OIL POLLUTION CAN BE DIVIDED INTO TWO GROUPS, 
MECHANICAL DAMAGE AND TOXIC EFFECTS. 



A, MECHANICAL DAMAGE. 



A MOST OBVIOUS EFFECT OF AN OIL SPILL IS THE AROUSED 
CITIZEN INTEREST WHEN NEARBY RESOURCES ARE POLLUTED BY THE 
SPILL. THE VISIBLE RESULT MAY BE AN OILED BEACH WHICH HAS 
BECOME LESS ATTRACT I VE , C I SCOURAGI NG SWIMMING AND FISHING. 
ALTHOUGH OIL COATED BEACHES CAN BE CLEANED COSMETICALLY, 

THEY CANNOT BE DISINFECTED (REF-2). 

IF LARGE CUANTITIES ARE SP I LL ED , TH ES E EFFECTS CAN BE IMPORT- 
ANT. SEA BIRDS ARE PERHAPS THE MOST OBVIOUS SUFFERERS. 

THEY DIE WHEN OIL DESTROYS THE NATURAL INSULATING QUALITIES 
OF THEIR FEATHERS (REF-8) OR MAKES FLIGHT IMPOSSIBLE, RESULT 
ING IN DEATH BY STARVATION. DUCKS WERE FOUND TO PREEN ABOUT 
HALF THE POLLUTING OIL FROM THEIR FEATHERS WITHIN A WEEK, 
MOSTLY ON THE FIRST DAY OF OILING . LATER EXPERIMENTS 
( HARTUNG R. AND HUNT, G.S.J WILDE. MGMT, 1966, P564-570) 
SHOWED THAT THE INGESTION OF ONE THIRD TO HALF THE AVERAGE 
QUANTITY OF OIL EXTRACTED FROM THE PLUMAGE OF BIRDS FOUND 
DEAD FROM MODERATE OILING, PRODUCES SERIOUS INTERNAL EFFECT- 
S. BUT THE IMMEDIATE CAUSE OF DEATH IN A BIRO OILED AT SEA 
IS 5XPCSLRE OR DROWNING. 

IN ADDITION, OIL DISRUPTS THE BIRDS' NESTING GROUNDS. 

ASIDE FROM BIRDS, SIGNIFICANT IMPACT DOES NOT GENERALLY 
OCCUR UNTIL OIL REACHES THE INTERTIDAL ZONE. 

THE MOST ecologically IMPORTANT IMPACT COULD OCCUR ON THE 
SHORE OF MARSHES AND ESTUARIES WHICH HAVE EXTREMELY HIGH 
BIOLOGICAL PRODUCTIVITY (REF-1,2,7). 

SMOTHERING OF SHORE CRABS AND SEiSILE INVERTI BRAT ES (NON 
COMMERCIAL WATER ANIMALS) CAN OCCUR IN THE UPPER INTERTIDAL 
ZONE OF THE OPEN COAST, BUT ARE OF CONSIDERABLY LESS ECOLO- 
GICAL SIGNIFICANCE (BRUCE SEYAERT, ANALYSIS OF OIL ACCIDENTS 
FOR ENVIRONMENTAL IMPACT STATEMENT, REF-1 P 43). 

ON THE SHORE ITSELF MOST SEAWEEDS HAVE AN OUTER LAYER TO 
WHICH OIL DOES NOT CLING, SO MODERATE QUANTITIES CF OIL CAN- 
NOT DC MUCH CAMAGE. INDIRECTLY IT CAUSES DAMAGE TO OTHER 
SPECIES WHICH DEPEND ON THESE PLANTS AS FOOD. 

A HEAVY SPILL WILL BLANKET EVERYTHING ON THE SHORE AND CLOG 
THE GRASSES AND REEDS OF S ALT-M ARSHcS , P ART I CUL ARL Y IF THE 
OIL HAS BECOME EMULSIFIED. A HEAVY RESIDUE CF CRUDE OIL 
REDUCED THE POPULATION OF WINKLES (LARGE SEA SNAILS) 

ALTHOUGH THE RESIDUE WAS VIRTUALLY NON TOXIC. 

A COATING OF OIL, EITHER ON THE SURFACE OF WAT=R OR CN AN 
INDIVIDUAL PLANT, ALSO INTERFERES WITH LIGHT PENETRATION AND 
THUS PHOTO-SYNTHESIS (REF-7). 



B. TOXIC EFFECTS. 



THE PHYTOTCXIC EFFECT OF HYDROCARBONS HAD BEEN STUDIED 
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IN MOST DETAIL ON TERRESTRIAL PLANTS . IT WAS FOUND THAT 
TOXICITY IN SMALLER MOLECULES WAS GREATER THAN IN LARGER 
ONES. 

FISH USUALLY KEEP WELL CLEAR OF AN OIL SPILL IF THEY CAN. 
SHELLFISH ARE AFFECTED BY OIL DUE TO THEIR LACK OF MOBILITY. 
FIVE TO TEN PERCENT CRUDE OIL SLOWS THE PUMPING RATE OF 
OYSTERS THUS AFFECTING THEIR FEEDING, RESPIRATION AND GENE- 
RAL CONDITION (REF-7). 

THIS HAS MEANT FINANCIAL LOSSES FOR FISHERMEN AND PROCESSORS 
. ALSC, SOME COMMERCIAL SPECIES CAN ACCUMULATE POTENTIALLY 
CARCENOGENIC SUBSTANCES, DAMAGING THE ORGANISM ITSELF OR 
MAKING IT UNFIT FOR CONSUMPTION BY MAN AND OTHER ANIMALS. 
AQUATIC LIFE UNDER CONDITIONS OF LONG TERM OR CONTINUOUS 
EXPOSURE TO OIL SPILLS OEVcLOPES SUBTLE CHANGES IN THE BEHA- 
VIOR PATTERN SUCH AS LOSING THEIR ABILITY TO SECURE FOOD, 
AVOID INJURY, ESCAPE FROM ENEMIES, CHOOSE A HABITAT, RECOG- 
NIZE TERRITORY, MIGRATE, COMMUNICATE AND REPRODUCE (REF-8) . 



C. IMPACT OF A DISCHARGE. 



AN OIL SPILL CAN HAVE ADVERSE SOC lA L, ECONOM IC AND EN- 
VIRONMENTAL IMPACT . THE SEVERITY OF THE IMPACT DEPENDS 
ON THE ENVIRONMENTAL SETTING OF THE AFFECTED AREA, THE TYPE 
AND AMOUNT OF OIL SPILT, AND THE MITIGATION MEASURES EMPLOY- 
ED. 

SOCIO-ECONOMIC IMPACTS COULD INCLUDE ADVERSE PUBLICITY AND 
PUBLIC CCNCERN WHICH COULD TEMPORARILY DETER COASTAL RECREA- 
TION AND TOURISM. 

PERSONAL PROPERTIES SUCH AS BOATS AND FISHING GEAR COULD BE 
FOULED BY OIL. 

COMMERCIAL AND SPORT FISHING ACTIVITY COULD BE TEMPORARILY 
STOPPED OR CURTAILED IN THE AREA. 

THE POSSIBLE IMPACTS OF A DISCHARGE CAN INCLUDE (R5F-2) : 



A. HAZARD TO HUMANS THROUGH EATING CONTAMINATED SEAFOOD, 

B. DECREASE OF FISHERY RESOURCES, AND DAMAGE TC WILDLIFE 
SUCH AS SEA BIRDS AND MARINE MAMMALS, 

C. DECREASE OF AESTHETIC VALUE DUE TO UNSIGHTLY SLICKS 
OR OILED BEACH, DECREASING THE VALUE OF PRIVATE PRO- 
PERTIES AND RECREATIONAL ACTIVITIES AND TOURISM, 

0. DECREASE IN DIVERSITY AND PRODUCTIVITY OF SPECIES IN 
THE POLLUTED AREA. 

e. MODIFICATION OF H A BIT ATS , DEL AY I NG OR PREVENTING RE- 
CGLCNIZATION. 
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V. THE ECCNQi^ICS OF OIL SPILLS 



I\ M/iRCH 1967 THE OIL TANKER TGRREY CANYON FOUNDERED OFF 
THE SOUTHERN COAST OF ENGLAND, SPILLING 119,000 TCNS OF 
CRUDE OIL. THE OIL SLICK QUICKLY SPREAD ACROSS NEARBY WATERS 
ANC FOULED LARGE AREA OF ADJOINING ENGLISH AND FRENCH 
COASTS. THE BRITISH GOVERNMENT ALONE SPENT SB MILLION (REF- 
12) CN CLEAN UP. THAT WAS ONLY A PORTION OF TOTAL CLEAN UP 
COSTS. IN ADDITION, THERE WAS EXTENSIVE LOSS OF MARINE LIFE 
AND FOULING OF BEACHES AND COASTLINES. 



SINCE OIL IS NOT COMPLETELY BIODEGRADABLE OR DOES NOT DETER- 
IORATE RAPIDLY, SLICKS AND GLOBULES OF OIL ARE VISIBLE 
THRQLGHOUT THE HIGH SEA OF THE WORLD. THE EXACT BIOLOGICAL 
CONSEQUENCES ARE STILL UNDETERMINED. 



IT IS TE 
SHOULD 6 
COME IF 
THE BENE 
FROM THI 
INDEED 8 
IT IS CR 
GET E PM IN 
PRODLCTI 
OF OIL S 
SECTION 



MPTING FOR 
E STOPPED, 
THE COSTS 
FITS. 

S POINT OF 
£ SOCIALLY 
ITICAL TO 
E THE SOCI 
VE INQUIRE 
PILLAGES. 
ARE DEVOTE 



PEOPLE TO ASSERT THAT ALL POLLUTION 
BUT THE SOCIETY WILL HAVE LESS REAL IN- 
OF TOTAL ELIMINATION OF POLLUTION EXCEED 

VIEW SOME LEVEL OF OIL POLLUTION MAY 
DESIRABLE. 

DEVELOP A PROCEDURE AND METHODOLOGY TO 
AL COST OF AN OIL SPILL SO AS TC CONDUCT 
S IN DERIVING THE SOCIALLY OPTIMAL LEVEL 
THIS SECTION AS WELL AS THE FCLLOWING 
D TO SUCH DEVELOPMENT. 



IF A SPILL OCCURS, THE DIRECT LOSS OF A PRODUCT TO THE ECO- 
NOMY MAY BE MEASURED IN TERMS OF THE MARKET VALUE OF THE 
PRC CLCT . 

THE^INDIPECT COST ASSOCIATED WITH THE ENVIRONMENTAL DAMAGE 
IS MUCH MORE COMPLICATED. 

THE INCIFECT COST IS THE FUNCTION OF SIZE, FREQUENCY, 
LOCATION AND TYPE OF OIL SPILLED INTO THE WATERS. 

THE SOCIAL COST IS DEFINED AS VALUATION OF LOSSES IN REAL 
GOODS ANC SERVICES RESULTING FROM THE CIL SPILL. 



IN THE ABSENCE OF ANY CLEAN UP PROCEDURES (DETECTION, 
CHEMICALS EQUIPMENT ETC), THE SOCIAL CCST OF A SPILL COULD 
BE DEFINED AS THE SUM OF DIRECT AND INDIRECT COSTS. 



C(S) = A(S) + B(S) 



WHERE: 

A, THE DIRECT COST IS A FUNCTION OF SPILL SIZE AND 

e, THE INDIRECT COST IS THE FUNCTION OF SIZE, 
FREQUENCY , LOCATION AND TYPE CF OIL. 

TYPE CF OIL CHARACTERIZED THE TOXICITY TO MARINE 
LIFE. 



A. GENERAL ECONOMIC ANALYSIS 
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ASSUME that there ARE N COMMOOTTIES IN THE ECONOMY WHICH 
CAN SERVE AS GGGOS WHERE THE COMMOOITItS APE DEFINED FOR A 
PARTICULAR DATE AND PLACE 30 THAT A SINGLE PHYSICAL COMMO- 
DITY OELIVEREC AT TWO DIFFERENT DATES CR TWO DIFFERENT 
PLACES WOULD BE CONSIDERED DIFFERENT ECONOMIC CCM '^GOI T I E S . 
ASSUME N IS FINITE, AND THE QUANTITIES OF ANY CGMVOCirY ARE 
assumed PERFECTLY DIVISI3LE. 

A PARTICULAR BUNDLE OF COMMODITIES IS SUMMARIZED BY THE 
COLUMN VECTOR *X': 



X = ( Xj^,X2»...»X^ )• 

THIS VECTOR IS DEFINED ON EUCLIDEAN N-SPACE,- i REFEREC 
TO AS CCMMCDITY SPACE. 

PRICE IN THE ECONOMY ARE SUMMARIZED BY ROW VECTOR P : 



P — ( P f ^ 2 ^ 



PRICE are non (NEGATIVE AND AT LEAST ONE PRICE IS NON ZERO. 
THE PRICE CAN BE NORMALIZED AND ONE POSSIBLE N OP M AL I Z AT I ON 
IS THAT OF MEASURING PRICES SO THAT THEY SUM TO UNITY. 



N 



EACH CF THE FIRMS IN THE ECONOMY MUST SELECT LEVELS OF 
INPUTS AND OUTPUTS SUBJECT TO THE AVAILABLE TECHNCLCGY, SO 
AS TG MAXIMIZE PROFITS. FOR EXAMPLE THE FIRM f MAY CHOOSE 

INPUT-OUTPUT VECTOR Y^ IN THE COMMODITY SPACE : 



THE FRGCLCTICM POSSIBILITIES SET A SUBSET CF COMMODITY 

SPACE IN which: 






f =1,^ 



IT IS ASSUMED THAT £aCH PRODUCTION POSS I SIL I T I ES SET IS IN- 
DEPENDENT OF THE INPUT/OUTPUT VECTOR CHOSEN BY CTHER FIRMS 
AND OF THE CONSUMPTION CHOICES OF CONSUMERS. 

THE ECONCMY -WIDE I-O VECTOR Y , IS OBTAINED BY SUMMING 
ALL INDIVIDUAL FIRM I-O VECTORS: 



Y 




F f 
SY 

{•I ^ 




) 



BY SUMMATION, INTERMEDIATE GOODS CANCEL OUT, SO ONLY FINAL 
QUTPLTS (MEASURED AS POSITIVE) AND PRIMARY RESOURCES (MEA- 
SURED AS NEGATIVE) APPEAR IN Y , 

THE ECONCMY-WIOE PRODUCTION POSSIBILITIES SET T IS OBTAINED 
BY SLMMING ALL FIRM PRODUCTION PGSS I B I L IT I ES SETS : 



F_ £ 

Y € T = IT 

f«t 
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ASSU^PTICNS : 

1. T IS CONVEX: Y,Z > ctY + (1- a)Z et 0 < a< 1 

2, IT IS IMPOSSIBLE TO PRODUCE OUTPUTS USING NC INPUTS. 

2. OUTPUT AND INPUT CANNOT BE REVERSED. 

4. IT IS POSSIBLE TO USE INPUTS ANO PRODUCE NO OUTPUT. 
INPUTS BEING FREELY DISPOSABLE. 

SINCE OUTPUTS ARE MEASURED AS POSITIVE AND INPUTS AS 
NEGATIVE, THE PROFIT OF FIRM f : 

f - -f f 

n = p.Y^ = 2P. .Y; 

J-t J J 



TOTAL PROFIT n IS MAXIMIZED WITHIN f IFF ALL FIRMS MAXIMIZE 
THEIR INDIVIDUAL PROFITS WITHIN THEIR PRODUCTION POSSI- 

BILITIES SETS . 

EACH OF THE CONSUMER H IN THE ECONOMY MUST SELECT LEVELS 
OF PURCHASES SUBJECT TO A BUDGET CONSTRAINT. 

CONSUMER h SELECTS A CONSUMPTION VECTOR C^ : 



THE TASTE OF CONSUMERS ARE SUMMARIZED BY THE PREFERENCE 
RELATION, ASSUMED CONVEX ANO CONTINUOUS. ALSO ASSUMED THAT 
THE PREFERENCE RELATION FOR ANY CCNSUMERS IS INDEPENDENT OF 
THE CCNSUMPTICN CHOICES OF CTHER CONSUMERS. 

THE EUCGET CONSTRAINT 3 : 



- V. ^ u 

B = P.C = 2P .C^ 
J=i j j 



total consumption levels for the ECONOMY C ,IS OBTAINED 8Y 
SUMMING ALL INDIVIDUAL CONSUMER CONSUMPTION VECTORS : 



Kii 



IF A REPRESENT TOTAL RESOURCES FOR THE ECONOMY, THEN 
WEALTH OF THE ECONOMY W : 



W = P.A 



A COMPETITIVE EQUILIBRIUM IS DEFINED AS A SITUATION IN WHICH 
PRICE VECTOR SATISFIES : 
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L 



<t.(p*) < 0 



i = 1,2, N 



(non positive excess demand function) 



P'l * 4>i ( P*) = 0 



t = 1,2,. . . ,N 



THE PROFIT MAXIMIZING I -0 VECT0=^ OF EACH FIRM IS SUMMARIZED 

_* _1* ' _2* _F* 

Y =: ( Y ,Y , . . . ) 

THE EQUILIBRIUM CONSUMPTION VECTOR OF EACH CONSUMER : 

_* — H* 

C = ( C ,C , . . . ,C ) 

PR’CFIT maximizing subject to THE AVAILABLE TECHNOLOGY AND 
POLLLTICN LEVEL : 



_* _f* _* _f 

= P .Y > P ,Y 



FOR ALL Y^ € T 



A PARETO .GPTi MUM IS A SET OF CONSUMPTICN VECTORS : 



1 * 2 * 

( C » C » 




) 



WHICH. IS CONSISTENT WITH THE TECHNOLOGY AND BUDGET 
AND FOR WHICH THERE EXISTS NO OTHER SET CF CONSUMPTION 
VECTORS, SUCH THAT NO CONSUMER IS WORSE OFF AND AT lEAST 
ONE IS SETTER OFF. 

IN THE CASE OF ONE CONSUMER AND ONE PRODUCER IS ILLUSTRATED 
IN FIG- 7. 







FIG-7 
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BY ASSU^PTICM THE PREFERENCE SETS AND THE PROOUCTIGN POSSIEI 
LITIES ARE CONVEX, COMPETITIVE EQUILIBRIUM IS GIVEN BY THE 
POINT OF TANGENCY OF THE BOUNDARY OF THE PRODUCTICN FP^ONTIER 
AND THE HIGHEST ATTAINABLE INDIFFERENCE CURVE, WHERE THE 

-h 

VECTOR GF THE CONSUMER IS C AND THE I-O VECTOR OF THE FIRM 

_f 

IS Y • 

BY CONVEXITY ASSUMPTIONS THERE EXISTS A SEPARATING HYPER- 
PLANE FOR WHICH THE PRODUCTION POSSIBILITIES SET LIES ON 
ONE SIDE AND THE PREFERENCE SET ASSOCIATED WITH THE HIGHEST 
ATTAINABLE WELFARE CURVE LIES ON THE OTHER SIDE OF THE 
HYPERPLANE. THE HYPERPLANE OR PRICE LINE FUNCTION : 



P .Z = V 
WH ER E J 

_• — * 

V = P .C = P .Y 

IN OUR STUDY WE ASSUMED THE X VECTOR COMMODITIES REFERS TO 
GOODS AND SERVICES SENSITIVE TO THt OIL POLLUTION SUCH AS 
FISHING INDUSTRIES, RECREATION ACTIVITIES ETC. 

ASSUMING THAT THE PRODUCTION POSSIBILITIES SET IS A FUNCTION 
OF SPILL SIZE. WE THEN ASSUME INCREASES IN POLLUTION SIZE 
WOULD CECPEASc THE PRODUCTION OF COMMODITIES IN THE mREA 
FOULED BY OIL. THIS MAKES SENSE SINCE FOULED BEACHES MEAN 
DECREASES IN INPUT GF RECREATION PRODUCTION, AND FOULED 
WATERS RESULT IN REDUCED COMMERCIAL FISHING REVENUES. 

IF GOODS AND SERVICES AFFECTED BY THE OIL SLICK HAVE NO 
MARKET VALUES OR PRICES, W5 MAY USE A SHADOW PRICING SCHEME 
TO DETERMINE THE VALUES. 

THE SHADOW PRICE OF A GIVEN (CONSTRAINED) COMMODITY IS DEFI- 
NED AS MARGINAL VALUATION OF THE COMMODITY BY RELAXING THE 
CONSTRAINT MEASURED IN TERMS OF THE OBJECTIVE FUNCTION . 

THE SOCIAL COST C(S) OF THE OIL POLLUTION COULD THUS BE 
DESCRIBED AS THE DIFFERENCE IN VALUE BETWEEN THE PRODUCTION 
LEVEL IN THE ABSENCE OF SPILLS AND THE PRODUCTICN LEVEL 
RESULTING FROM A SPILL. 



C(S) = P(S=0) 



x1s=0) - X*S = 1) ) 
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VI. COST MODEL DEVELOPMENT 



ANY SPILLAGE OF OIL INTO THE WATERS REPRESENTS A LOSS OF 
PROCLCTS TO THE ECONOMY. THIS LOSS IS EQUAL TO THE MARKET 
VALUE CF THE PRODUCTS LOST, 

IN ADDITION TO THE LOST PRODUCT, THERE WILL ALWAYS BE DAMAGE 
DONE TO THE AQUATIC ENVIRONMENT AND POSSIBLY TG THE SURROUN- 
DING LAND AREAS OR BEACHES. 

MANY FRACTIONS OF OIL ARE SOLUBLE OR EMULSIFIED IN THE WATER 
AND SINK. SUCH DISSOLVED SUBSTANCES CANNOT BE EASILY TAKEN 
OUT EY SIMPLE REMOVAL OF THE VISIBLE PRODUCT AND DAMAGES 
INFLICTED MAY NOT BE REVERSIBLE. 

THE ENVIRONMENTAL DAMAGE FACTOR IS COMPLICATED BY THE 
EXISTENCE OF SHORT AND LONGER TIME EFFECTS. 

THE LOCATION OF THE SPILL IS CRITICAL TO THE COST THAT IS 
INCURRED. 



A. GENERAL COST MODEL. 



AS WAS PREVIOUSLY MENTIONED, THE COST OF A PARTICULAR GIL 
SPILL C(S) : 



C ( S ) = A( S ) 




WHERE: 

A (S) = DIRECT COST OF OIL SPILL. 

B (S) = INDIRECT COST OF OIL SPILL. 

8. DIRECT COST, 

DIRECT COST (S) IS EQUAL TO THE MARKET VALUE OF THE 
PRODUCT LOSS 

A (S) = C .V 



V = VOLUME OF GIL SPILLED 
C = MARKET PRICE PER UNIT VOLUME 
S = SPILL INDICATOR. 



C. INDIRECT COST. 
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INDIRECT COST B(S), IS THE VALUATION OF COMMOCITY LOSS 
RESULTING FROM THE OIL SPILL. 



N 



B ( S) = C. Ax = ^ C; . 






TO SIMPLIFY THE COST MODEL TO A RESOLVABLE LEVEL, CONSIDER 
ONLY COMMODITIES THAT MIGHT HAVE SIGNIFICANT EFFECT IN THE 
CALCULATION OF COST. 

ASSUME THE SIGNIFICANT INDIRECT COST WAS THE SUMMATION OF 
LOSS IN WILD BIRD POPULATION, LOSS IN FISH POPULATION, LOSS 
IN PROPERTY VALUES, LOSS TO THE BUSINESS ACTIVITIES IN 
COMMERCIAL FISHING AND RELATED INDUSTRIES AND LOSS TO THE 
RECREATION ACTIVITIES ( SPORT-FISHING, BOATING, TOURISM ETC) 



THUS IN CUR MODEL COST, N EQUAL TO FIVE, WE DENOTE : 



= RELATIVE PRICE (WEIGHTING FACTOR) CF A BIRO. 
AXi = NUMBER OF BIROS KILLED. 



THE RELATIVE PRICE REFLECTS THE SHADOW PRICE WHICH IS PRE- 
SUMABLY DETERMINED UNDER SOME APPROPRIATE CONSTRAINED MAXIMI 
ZATICN PROBLEM. 



C2 = RELATIVE PRICE OF A FISH 
AX2 = NUMBER OF FISH KILLED 



AS IN C -1 ,USE SHADOW PRICE TO DETERMINE THE RELATIVE PRICE 
C-2 OF A FISH. IN THIS CASS THE FISH ARE NOT VALUED AT 
THEIR MARKET PRICE TO AVOID DOUBLE COUNTING, SINCE THE 
EFFECT ON THE MARKET IS COMPUTED BELOW AS THE REDUCTION OF 
OUTPUT IN the COMMERCIAL FISHING ACTIVITY. THE SHADOW PRICE 
REPRESENTS THE LONG TERM EFFECTS AND ALSO REPRESENTING THE 
real obvious loss in the NATURAL RESOURCE BASE. 



= DROP IN MARKET PRICE OF PROPERTY VALUE PER 
^ UNIT AREA RELATIVE TO THE BASE PRICE. 

AX5 = AREA OF VALUABLE B EAC HES /COAST FOULED BY OIL. 

C -4 IS A FUNCTION OF LOCATION. IN OUR STUDY WE DIVIDE SAN 
FRANCISCO BAY AREA INTO FIVE DIFFERENT LOCATIONS. THIS CCULD 
BE REFINED INTO ON THE SPOT LOCATION AS PREDICTED BY THE 
SIMULATION TECHNIQUE. IN REALITY WE SHOULD ALWAYS USE MUCH 
FINER GRID SYSTEM TO OBTAIN A MORE ACCURATE COST PICTURE. 



C4 = MARKET PRICE PER UNIT PRODUCT OF COMMERCIAL 
FISHING ACTIVITIES (AND RELATED INDUSTRIES) 
RELATIVE TO THE BASE PRICE. 

AX4 = REDUCTION ON PRODUCTION CAPACITY OF COMMER- 
CIAL FISHING ACTIVITIES (AND RELATED 
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INDUSTRIESJ 



C 



5 






MARKET PRICE PER UNIT PRODUCT OF RECREATION 
ACTIVITY RELATIVE TO THE BASE PRICE. 

REDUCTION IN PRODUCTION CAPACITY OF RECREA- 
TION ACTIVITY. 



•D. DETERMINISTIC MODEL TO ESTIMATE BIRD OR FISH DEATH. 



TO ESTIMATE THE BIRO OR FISH KILLED, WE USE DIFFERENTIAL 
EQUATION AS FOLLOWS: 



dX 

dt 



=. -k.Y 



WHERE : 

k = DEATH RATE OF BIROS OR FISH AS A FUNCTION OF 
LOCATION, SIZE AND TYPE OF CIL SPILL. 

TYPE OF OIL RELATES WITH ITS TOXICITY. 

Y = AREA SPILLED AS A FUNCTION OF TIME. 

SOLVE THE EQUATION BY USING EULER APPROXIMATION TECHNIQUE: 

NUMBER OF BIRDS OR FISH KILLED = AX 

M 

Axl 5 k.At. S Y. 

lt=T ^ 



E. DETERMINISTIC MODEL TO ESTIMATE PROPERTY DAMAGE. 



IF THE OIL SLICK LANDED ON THE SHORE, THUS THE SOCIAL 
COST OF DAMAGE WOULD BE CONSIDERED TO BE EQUAL TC THE LOSS 
OF VALUE IN MARKET PRICE OF THE PROPERTY FOULED BY GIL. 



2 

StTh^-. 

• A Xj 's C- .TT. 



2 

TTPf REPRESENT THE AREA FOULED BY OIL. Rf MAY BE OBTAINED 
FROM THE SIMULATION. 
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F. MODEL TO ESTIMATE THE SOCIAL COST OF BUSINESS ACTIVITY 



THE STANDARD PRODUCTION FUNCTION MODEL IS A SUFFICIENTLY 
WELL KNOWN TOOL OF ECONOMIC ANALYSIS. 

TO MAKE THE MODEL TRACTABLE, WE SIMPLIFY BY ASSUMING THAT 
ONLY A SINGLE FIRM POSSESSES ALL BUSINESS ACTIVITIES IN TH 
BAY AREA, AND THE EXISTENCE OF SOCIAL PREFERENCE OVER RECRE 
TION ANC COMMERCIAL FISHING ACTIVITIES. 

ASSUME THE PRODUCTION POSSIBILITIES SET OF THE FIRM IS 
KNOWN AS FOLLOWS : 




WHERE : 

X = PRODUCT OF RECREATION ACTIVITY 

Y = PRODUCT OF COMMERCIAL FISHING (AND RELATED IN- 
DUSTRIES) 

a ^ 0 , b ^0 

a ,b IS A FUNCTION OF SPILL SIZE AND LOCATION. 



c 2 

r ^ ■ f 

a = -H 6 ^ } 

-h 'F 

b = bo(^e 



(^1 = THE PROPORTION OF PRODUCTION RELATED WITH THE 

^ BAY ACTIVITIES ( SWIMMING, SPORT FISHING ETC) 

S 2 = DEGRADATION CONSTANT 

= THE PROPORTION OF PRODUCTION NOT RELATED TO THE 
^ BAY ACTIVITIES FOR EXAMPLE PEOPLE IN TRANSIT ETC 

= THE PROPORTION OF PRODUCTION WITH LOCAL INPUT 
RESOURCES. 

= DEGRADATION CONSTANT 

fX = THE PROPORTION OF PRODUCTION WITH INTERLOCAL 
5 INPUT RESOURCES. 

V = VOLUME, SPILL SIZE. 

8U , = MAXIMUM PRODUCTION LEVEL FOR EACH ACTIVITY 

° GIVEN ALL AVAILABLE RESOURCES. 

a^, bp, 5,6 IS A FUNCTION OF LOCATION. 
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l>rn 



THE SOCIAL PREFERENCE FUNCTION IS GIVEN AS 



W(s; = c.x'^.Y^ 



WHERE : 

C = CONSTANT, GIVEN AS A FUNCTION OF LGCATIGN. 

C<, : KNOWN AS GIVEN AN0CK+P<1»<^>0, (3>0. 

BY ASSUMING THAT THE FIRM TREATS PRICES AS GIVEN ANC SEEKS 
TO MAXIMIZE ITS PROFIT, 

MAXIMIZE ITS PROFIT, 



MAX = P .X + P^ .Y 
X A' 



SUBJECT TO 




0 



BY USING LAGRANGE MULTIPLIER TECHNIQUE, THEN 



PROFIT IS MAXIMIZED IF ; 




\ 



ON THE CCNSUMER SIDE, GIVEN BUDGET CONSTRAINT, MAXIMIZE HIS 
PREFERENCES 

MAX VV(S) = C.X .Y 

SUBJECT TO Z = P .X + P ,Y 

USING LAGRANGE MULTIPLIER TECHNIQUE : 
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PREFERENCeS ARE MAXIMIZED IF 



o< -z. 



EQUILIBRIUM CONDITIQN IF: Z = TC = Tt 

THE CONDITION APPLIES IF : 





THE SUPPLY ANC DEMAND CURVES OF 'X' PRODUCT ARE 



a'-p* 

x, = -^ 

^ TT 
V c<rr* 

■ 



THE SUPPLY AND DEMAND CURVES OF 'Y* ARE : 
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2 * 
b^Py 

TT * 



Yc 



IS.TT' 



Py(^+ p) 



WHERE : TL = 



/ 2^*2 , ^2 *2 , 1/2 
( a Pjj ^ ) 



THE PRICE FUNCTION IS : 



= P 



Y . a 

e 

2 2 1/2 
b( b - ) 



G. SENSITIVITY ANALYSIS. 

AS AN EXAMPLE WE ANALYSE THE CHANGES IN SUPPLY CURVES 
AS WE CHANGE THE PRODUCTION POSSIBILITY PARAMETER : 

CHANGES IN SUPPLY CURVES : 



■aa 

Tia 



a.b^ (P /P 
X y 



a 



a 



a 



( 



2 2 2 
(P^/Py)^ + b^ ) 

^ + 2b^ (Py/?x ) 

2 2 2 
a + b (Py/P^ ) 



3/2 

2 



4.0 




3 5 



CONCLUSICN 



1. IF 'A* DECREASES, THE SUPPLY CURVE OF WILL 

SHIFT TO THE LEFT, IF *A' INCREASES CONVERSE HOLDS 

2. IF ‘A* DECREASES THE SUPPLY CURVE OF 'Y* WILL 
SHIFT TO THE RIGHT, IF 'A* INCREASES CONVERSE 
HOLDS. 



CHANGES IN DEMAND CURVES : 



CONCLUSION : 

1. IF 'A‘ DECREASES, BOTH DEMAND FOR *X‘ AND 'Y* WILL 
SHIFT TO THE LEFT. 

2. IF 'A* INCREASES, CONVERSE HOLDS. 

CHANGE IN THE EQUILIBRIUM POINT : 






2>a 







■36 /-^ . a • 




('<=!'+ |3) 
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=. o 



^ Xc. _ 2. <3 oC ^ 3 ^ 0 ^ \ '■^2. 

'2>a ■ o< + p, ^ c?<4- 
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CQNCLUSICN 



1. NO CHANGE IN *Y' PRODUCTION LEVEL. 

2. IF 'A' DECREASES, THE PRODUCTION LEVEL OF ‘X* WILL 
DECREASE, IF »A* INCREASES , CONVERSE HCLDS. 

THE SOCIAL COST OF A SPILL IS EQUAL TO THE MARKET VALUE OF 
THE FRGCUCT LOSSES. 



SOCIAL COST = P^.{Xg- X'g ) + P^ . ( - Y^' ) 



WHERE : 

Px AND Py REPRESENT THE PRICE PER UNIT PRODUCT OF X 
/IND Y, 

X„ AND Y« REPRESENT THE OPTIMUM PRODUCTION LEVEL 
OF X AND Y BEFORE THE SPILL OCCURS. 

x' AND y' represent THE OPTIMUM PRODUCTION LEVEL OF 
0 a 

X AND Y AFTER THE SPILL OCCURS. 



SEE FIGURE 8A,B ; FIGURE 9A, ; FIGURE 10A,B FCR GRAPHICAL 

ILLUSTRATION. 



37 




X UNIT PRODUCED 



- PROD. ^033 CURUe 

- hPIK PROFIT LINE 
* UELFRRE CuRUE 





\ \ * 1 


fig-8b 





a * 

K UNIT PRODUCED 




K UNIT PRODUCED 




K UNIT PRODUCED 




V UNIT PRODUCED 




y UNil PROCUClO 



VII 



DATA 



ANALYSIS 



ASSUME THE PARAMETEi^S OF THE COST MODEL IN THIS STUDY ARE 
SHOWN IN TABLE-2 



REGION 


PROD. PARM 


BEACH VA 


,LU£ 


BIRD PARM 


FISH PARM 


SOUTHERN 


=40,000 


H=4,000 


N 


C^= 2 


N 


€ 2 = 2 N 


BAY / EAST 


=10,000 


L=2,500 


N 


k^=o:oi 




k2 =0.12 


SOUTHERN 


aQ =60,000 


H=6,000 


N 


Ci= 3 


N 


C2= 2 N 


BAY / WEST 


bg =40,000 


L=3,500 


N 


lcj^=0.01 




k2 =0.10 


CENTRAL 


a^ =100,000 


H=7,000 


N 


Cl= 5 


N 


C2 = 2 N 


BAY 


bg = 50,000 


L=3,500 


N 


ki=0.05 




k2 =0.05 


NORTHERN 


=20,000 


H=4,000 


N 


Ci= 3 


N 


{^ = 4 N 


BAY / EAST 


bo =30,000 


L=l,500 


N 


0 

• 

0 

It 

1 — 1 
.M 




k2 = 0 . 1 5 


NORTHERN 


aQ =100,000 


H=8,000 


N 


c^= 10 


N 


(^ = 1C N 


BAY / WEST 


bo =100,000 


L=5,000 


N 


k^ =0.10 




k2 =0.15 


* N=NUMERAIRE ; BASE 


PRICE, ONE 


LIT 


ER CIL = 


ONE 


NUMERAIRE 



=!« 6EACF VALUE : N PER HECTARE. 

* BIRD (FISH) VALUE : N PER UNIT KILLED. 

P4RM = PARAMETER 

* H a L : H STANDS FOR THE VALUE OF BEACH OR OTHER 

property PRIOR TO THE OIL SPILLAGE AND L 
AFTER THE SPILL, THUS THE MARKET VALUE OF 
THE LOSS (C-5 ) IS GIVEN BY : 

C5 = H - L. 



TABLE-2 
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1 



THE BECPEiTICN S COMHSRCIIL FISHING COST MODEL PASArST'SRS; 



i, 


= O.S 


S., = 0.75 


4 


II 

O 

\ 


= 25X10"'^ 




= 0*1 


£3 = 0.25 


di 


= 0.45 


II 

0 

• 

u» 


Px 


= 150. N 


p = VARIES AS A FUNCTION CF 
V 


COG 


= COCROINAT 

latitude 

longitude 


c, COMPUTED IN METERS FROM THE ORIGIN 
" ; 37» 25' N 

: 122^ 35’ W 



USING T^-eSE fpIf3Sj^.6,g3TA!NED^F|OH 

INDUSTRIES) ACTIVITIti TAeULAi-D 



TIGN 

(AND 



VALUES 

related 



IN TAELE-5« 



LOSS IN RECREATION 


ACT. AND 


COMMERCIAL FISHING 


ACT. 








VOLUME (LITE 


RS) 


LOCATION 




lOtOOO 


100,000 


1 9 00 C t 0 0 c 


SOUTHERN 


BAY/EAST 


6,671. 


595,255. 


5,099,228. 


SOUTHERN 


bay/west 


10 ,005. 


392,382. 


7,648,843. 


CENTRAL 


EAY 


16,67 1. 


1,488,135. 


12,746,06C. 


northern 


BAY/V^EST 


16 ,673. 


1 t433t 133* 


12,748,066. 


NORTHERN 


BAY/EAST 


3,335. 


297,626. 


2,549,613. 






TABLE 


- 3 





TS 

ACCI- 

WILL 



TC DETAIN THE SOCIAL COST OF -k kV - t pm 

Ff irTli I? iffuSr/SE'S^r?YblL^S^fi?’ SNi^fHl^J^FLlTI 

SIX =xil-PLE^CF°COMPUTING THE SOCIAL COST WAS OETAINEO AS 
FQLLCWS : 



GV£ VEN' 
ON OF 
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SAMPLES OF COMPUTING THE SOCIAL COST USING THE COMPUTER 
SIMULATION TECHNIQUE: 



1 . COC : 


(20 200 *,50700 ) 




^.2 


= 0.4 


HRS 


VOL : 


100,000. LTRS 






= 1.9 


HRS 


sj = 


398 METERS 






= 5.8 


HRS 


LQCATICN 


AREA FLD 


BEACH FLD 


DEAD 


BIRD 


C5A0 FISH 


CENTRAL EAY 4234.1 


49.3 


212 




635 



SOCIAL COSTS = 1,764,765. N 



2. CCO : (20200;50700) T, 2 . = 0.9 HRS 

VOL : L,OOC,OCO. LTRS T 2.3 = S.9 HRS 

: 944 METERS = 13.3 HRS 



LCCATICN 




AREA FLD 


BEACH 


FLO 


DEAD 


BIRD 


DEAD 


FISH 


CENTRAL EAY 


6940. 


230. 




347 




1041 






SOCIAL COSTS = 


14, 731, 


877. 


N 








3. COO : 


(29000;64500 ) 








= 0.2 


HRS 




VOL : 


10 , 


COO LTRS 






Tz3 


= 0.4 


HRS 




Rj : 


168 


METERS 








= lt 8 


HRS 




LOCATION 




AREA FLD 


BEACH 


FLD 


DEAD 


BIRD 


DEAD 


FISH 


NORTHERN 

ViEST 


EAY 


29.4 


NONE 




2 




C 






SOCIAL COSTS = 


26 


,753. 


N 









^6 



, coo : (29C00;64500 ) 
VOL ' ItOOOfOOO LTRS 
Rj : 94A METERS 



Tgj =8.9 HRS 
= 18.3 HRS 



location 



AREA FLO BEACH FLO CEAO BIRD CEAD FISH 



JORTFERN BAY 2089.6 
WEST 

SOCIAL COSTS = 

5. CCC : (36000;21000 ) 

VOL : IfOCOtOCO LTRS 
Rj : 944 METERS 



78.3 209 313 



13»988»186. N 

T,^ = 0.9 HRS 

T 23 =8.9 HRS 

Tj = 18.3 HRS 



LOCATION AREA FLD BEACH FLO 



SOUTHERN BAY 2»432. NONE 

WEST 



CENTRAL BAY llf705.5 280. 



DEAD BIRO 



24 



1,171 



CEAO FISH 



243 



11,706 



SOCIAL COSTS = 22,406,971. N 

6 . COO : (37500;1S500 ) 

VOL : 100,000 LTRS 
Rj : 598 METERS 



T = 0.4 HRS 
12 

T^3 = 1.9 HRS 
= 5.8 HRS 



LOCATION area FLO BEACH FLD 



SOUTHERN BAY 1,640. 
EAST 



CEAD BIRO CEAD FISH 



16 197 
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SOUTHERN BAY 



3,168 



50 



32 



380 



SOCIAL COSTS = 1,714,799. 



THE AGGREGATE LOSSES SHOW THE SIGNIFICANT DIFFERENCE BET- 
WEEN SHALL ANC LARGE SPILLS. 

THE MOST SIGNIFICANT ECONOMIC LOSSES ARE SUFFERED BY THE 
BUSINESS ACTIVITIES IN THIS ILLUSTRATION. 

NO ATTEMPT HAS BEEN MADE IN THIS STUDY TO FORMULATE THE LONG 
TERM SOCIAL COST OF OIL SPILLS. 

THE MAGNITUDE OF COMPUTED VALUES IN THIS ILLUSTRATION NATUR- 
RALLY DEPEND UPON THAT OF THE PARAMETERS, WHOSE VALUES ARE 
ARBITRARILY CHOSEN, BUT ILLUSTRATIVE OF POSSIBLE PEAL WORLD 
SITUATIONS. 
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VIII, CONCLUSICN 



Th£ PURPOSE OF THIS THESIS IS TO DETERMINE THE SOCIAL 
COSTS OF OIL SPILLS. 

THE REALIZATION THAT SUCH COSTS DEPEND ON THE ORIGINAL LOCA- 
TION ANC SIZE OF A SPILL AND THE EVENTUAL AREA AFFECTED 
NECESSITATES DEVELOPMENT OF A MODEL WHICH PREDICTS HOW A 
GIVEN OIL SPILL WILL SPREAD UNDER PLAUSIBLE CIRCUMSTANCES. 
FOR THIS WE HAVE SYNTHESIZED THE WORKS OF FAY (REF-6), SMITH 
(REF-4), SIVADIER AND MIKOLAJ (REF-1) , CONOMOS (REF-5) BY 

INCORPORATING WIND (RANDOMIZED DIRECTION AND VELOCITY), 
ESTUARINE NON TIDAL DRIFT AND EVAPORATION FACTOR. 

IN ORDER TO DETERMINE SOCIAL COSTS OF AN OIL S'^ILL WE ALSO 
NEED TO IDENTIFY THE COST COEFFICIENTS OF FACTORS DAMAGED 
SUCH AS WILDLIFE, RECREATIONAL ANC COMMERCIAL ACTIVITIES, 
ETC. EXACT DETERMINATION OF THESE COEFFICIENTS IS BEYOND THE 
SCOPE OF THIS PAPER. THUS IN THE QUANTITATIVE ANALYSIS 
GIVEN HERE FOR ILLUSTRATION, THEY ARE GIVEN PLAUSIBLE, BUT 
ARTIFICIAL VALUES. 
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APPEND! X 



A 



COMPUTER SIMULATION PROGRAM OF 
THE SPREAD AND MOVEMENT OF OIL SLICK 
IN SAN FRANCISCO BAY 



//HAN187A JOB ( 2696, 1446 , RL44 ) t ' MUDJ I ARD JO SMC 1374* 

// EXEC FORTCLGP, REGION. GO=180K 
//FORT.SYSIN CD * 

THESIS THESIS THESIS THESIS THESIS THESIS THESIS THESIS TH 
THE SPREAD AND MOVEMENT OF OIL SLICKS ON THE WATERS. 

OSL = OIL SLICK LEEWAY 

WIND = WIND SPEED MEASURED AT 10 METERS Ei-EVATIGN IN KNOTS 

AN ESTIMATE OF WIND DRIVEN SURFACE CURRENT MAY BE OBTAINED 
FROM THE FOLLOWING EQUATION PROPOSED BY THORADS (1914) : 

sew = SURFACE CURRENT DRIVEN BY THE WIND 
PHI = geographical LATITUDE 
sew = 0.0361 * SQRT(WIND) /SORT ( (SIN( phi ) ) 

IF( WIND .GT. 11.64 ) SC W=0 .01 26*WIN0/SQRT ( ( SIN ( PHI ) ) 



DIMENSION TIMES(80)t RADIUS(80) 

OIMENSI ON X( 80) , Y( 80 ) 

DIMENSION CCAA(900), 0DAA(900) 

DIMENSION 0SLV(80), SCWV(30), WINDV(80) 

DIMENSION DELT( 80),HTT( 80) 

DIMENSION CC( SO, 30) , DD( 30, 80) 

DIMENSION CCA( 6400), QOA ( 6400) 

EGLIVALENCE(CCA,CC), ( CCA,0D) 

DIMENSION XMV(80), YMV(80) 

DIMENSION XP(200), YP ( 200 ) , XO ( 200 ) , YQ ( 2 00 ) 

DIMENSION DZ(IOO), HZ(IOO), VXL(80) 

DIMENSION XQ(200 ) , YQ( 200) ,XR(200) ,YR(200 ) 

DIMENSION XS (200 ) , YS (200) , XT ( 200 ) , YT ( 200 ) 

DIMENSION XU(200) ,YU(200) ,XV(200) ,YV(200) 

DIMENSION AX (800 ) ,AY( 300) 

96 FORMAT( 6F5.1 ) 

97 FQPMAT( 2F10.1 ) 

98 FGRMAT( F10.3 ) 

99 FORMAT( 513 ) 

100 FCRMAT(/,10X,'M=‘ ,I3,3X,'N=' ,I3,3X,'KK=' ,I3,3X,'NK = ', 
1I3,3X,'KX=‘, 13, 3X, ' XMV( 1 ) = ' , F3 . 1 , 3 X, • YM V ( 1 ) = ' , F3 . 1 ) 

101 FQRMAT(/,10X, 'VOLUME OF OIL SPILLED IN L IT E RS = ' , F 1 2 . 1 ) 

102 FORMAT ( lOX, • RADIUS OF SLICK, M',3X,'AT TIMS, HRS', 

23X, 'THICK OF 0 1 L , CM ' , 3X , ' TH I CK OF GIL P ART! CL ES , CM ' , 
33X, 'VOL IN LTRS' ,3X, 'NO: ' ) 

103 FORMAT ( 15X,F8.2, IIX, F6. 1 , 10X,F7.3, 15X, F7.3 ,12X,F10.1, 
45X,I3 ) 

104 FORMAT (//, 18F8. 2 ) 

105 FCPMAT( 18F8.2 ) 

106 FORMAT (/, lOX, ' MAXIMUM RADIUS OF OIL SLICK IN METERS=', 
5F8.2 ) 

107 FORMAT!///) 

108 FORMAT! • 4= 

109 FORMAT!/, lOX, ' CRITICAL THICKNESS IN CM:',F7.3 ) 

110 FORMAT! /, lOX, ' THICKNESS OF MOL. DIFFUSION, CM=',F7.3) 

111 FORMAT!/, lOX,' TIME REQUIRED AS INFINITY, HRS:',F6.1 ) 

112 FORMAT! '!' ) 

113 FORMAT ! lOX, 'GAMA 1= ' ,F5.1,3X,'GAMA2=',F5.1,3X,'GAMA3 = ', 
7F5. 1,3X, 'GAMA4=' ,F5.1,3X,'GAMA5=',F5.1,3X,' GAMA6=' , 
8F5.1 ) 

114 FORMAT! /, lOX, 'T12 = ',F6.1,3X,' T23 = ' , F6. 1 ) 
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2 01 FORMAT { lOX , F9 . 3 , 5X , F 9.3 , 3X , F 9. 3 , 2X , F9 . 2 , 4X , F9 . 2 , 6X , I 3) 
20 2 FORMAT ( 13X , » WI NO ' t 3X, • QSL • , 9X, • , 1 OX , » XM V • , iOX t 

1 ‘YMV* ,8X, ‘NO:* ) 

203 FORMAT ( / ) 

300 FORMAT ( 16F5.1 ) 

301 FGRMAT( 16F8.2 ) 

302 F0PMAT( 16F5.2 ) 

303 FGRMAT(/, lOX, ' RADIUS OF OIL =‘,F8.2 t* NZ= ',14 ) 

304 F0RMAT(/,10X, 'TOTAL AREA FOULED IN HECT A R E= * , F 14 . 1 ) 

READ 

GAMM4 = ANGLE BETWEEN THE WIND VECTORS AND THE COORDINATE 
SYSTEM 

READ(5,96) GAMAlt GAMA2t GAMAS t GAMA4, GAMA5, GAMA6 
SPILL COORDINATE MEASURED FROM ORIGIN IN METERS. 
REAO(5,97) XMV(l), YMV(l) 

GIVE THE INPUT DATA OF THE QUANTITY OF OIL SPILLED IN KILO 
LITERS. 

1 KILO LITER = 1000 LITERS. 

READ(5,98) QTT 

M,NtKK = FIXED NUMBER OF COORDINATE JOINTS CF SAM FRAN 
CISCO BAY AREA AND CURRENT VECTORS ON THE BAY AS PRECICTED 
BY ESTUARINE NONTIDAL DRIFT TO BE PLOTTED BY THE CGMOUTER. 
NK ; NUMBER OF ITERATION RELATED TO TIME NEEDED TO 
ESTIMATE THE OIL MOVEMENT. 

KX : NUMBER OF ITERATION TO COMPUTE THE COORDINATE OF 
THE OIL SPREAD. 

LIMIT OF (NK*KX) = 900 POINTS. 

READ(5t99) M, N, KK, NK , KX 
C COORDINATE POINTS TO CRAW SF BAY AREA. 

R5AD(5t300) ( XO (I ) , 1 = 1, KK ) 

READ(5,300) ( YO(I ), 1=1, KK ) 

C THE COORDINATE POINTS TO DRAW THE ESTUARINE NON TIDAL 
C DRIFT. 

READ(5,302) (AX(I),I=1,N ) 

REAOt 5, 302) (AY( I ) ,I = 1,N ) 

C the COORDINATE POINTS TO DRAW THE ESTUARINE NON TIDAL 
C DRIFT CONTINUED 

R6AD( 5, 302) ( XR ( I ) , I = 1 , M) 

REA0(5,302) ( YR ( I ) , I = 1 , M ) 

044444444444444444444444444444444444444444444444444444444444 

C 



C 

c 



Cl = 0.53 
C2 = 0.026 
ZU = 0.01 
RhCIL = 0.95 
RHOWAT = 1.0 
GRAV = 980.0 
D = 1.0 
SIGMA =10. 

PHI = 38. 

IX = 1234567 

OELRHG= RHOWAT - PHOIL 
DELTA= DELRHO/RHOIL 
VOL = GTT=^ 1000000.0 
VOLUME = VQL/1000. 



COMPUTE 



C COMPUTE THE CRITICAL THICKNESS OF THE OIL SLICK. 

HC = ( SIGMA/(OELRHO # GRAV))«-.5 
WRITE! 6, 112) 

WRITE(6, 108) 

WRITE(6,100) M, N, KK, NK , KX, XMV(l), YMV(l) 
WRITE(6,113) GAMAl, GAMA2, GAMA3, GAMA4, GA^AS, GAMA6 
WRIT£(6,101) VOLUME 
WRITE(6,109) HC 

C MAXIMUM RADIUS OF GIL SLICK AS ONE UNIT AT ZERO GROWS RATE 
C ( INDEPENDENT OT TIME ) 
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VVV = VQLUM=/1000. 

REND = 100. * ( 10.=!'*5*VVV*=!=0. 75/ 3. 14 )^*0.5 
CALL RANDU(IX,IY, YFL) 

C3 = 0.78 + YFL*0.02 
C ENC RADIUS REDUCED DUE TO EVAPORATION 
REND = RENO * C3 /2. 

TINF= ( REND/2.3 )=^*4*RHCWAT**2*ZU/S IGMA^^Z 

TINF = TINF*=i'0. 33333 

HINF = SQRT( TINF*D/10.*=S‘5 ) 

TINF = TINF/3600. 

REND = RENO/100. 

T12 = ( 1. 45/1.1 4 )**4*( VOL /(DELTA’f'GRAV^ZU) ) =!'*0.333 
T12 = T12/3600. 

T23=( 1.45/2. 30)*=«'2=!'(RH0WAT/SIGMA 1 * VOL=S'*0. 6667 
1( DELTA * GRAV ZU )**0.333 
T23 = T23/3600. 

VVL = VOL 
WRITE(6,106) REND 
WRITE(6,110) HINF 
WRITE{6tlll) TINF 
WRITE(6,114) T12,T23 
WRITE(6, 108) 

WRITE! 6,107) 

WRITE(6,102) 

TIMES! 1 ) ^ 0. 

TS = 0. 

DC 20 K = 1,NK 

C ASSUME THE WIND BLOWS RANDOMLY AT SPEED BETWEEN 0.5 TC 25 
C KNCTS PER HOUR. 

CALL RANDU!! X, lY, YFL) 

WIND = YFL 24.5 + 0.5 
XXX = SIN!0.01745*PHI ) 

C COMPUTE THE SEA SURFACE CURRENT VELOCITY INDUCED BY THE 
C WINDS. 

sew = 0.0361 * SQRT!WIND)/SQRT!XXX) 

IF! WIND .GT. 11.64 ) SCW = 0 . 0 1 26*WI ND/ SQR T ! XXX ) 

C COMPUTE THE OIL SLICK LEEWAY BLOWS BY THE WIND. 

OSL = 0.0199 =>= WIND 

IF!WIND.GT.5.0.AND.WIND.LT .25.0 ) OSL= . 0 1 79* wl NO+0 . 0 196 
C C3NVERT KNOTS/HR TO METERS/HR : 

WINDV!K)= 1850. * WIND 
SCWV!K) = 1850. *SCW 
OSLV!K) = 1350. *OSL 
C SPREAD BY GRA V IT Y- INE RTI A FORCE 
TIME = FLOAT!K) 

T = TIME*15,*60. 

TIME = TIME/4. 

IF (TIME .GE. T12) GO TC 1 
DELT!K) = !ZU*T)**.5 
VLL = VQL/10.**4 

R1 = 10.*! DELTA * GRAV * VLL * T**2 )**0.25*1.14 
HT = V0L/R1**2 
R1 = Rl/100. 

RADIUS!K) = R1 
TS = TS + 1. 

GO TO 21 

C SPREAD BY GRAVITY-VISCOUS FORCE 
1 AA = ZU**.5 

TIME = FLOAT!K) - TS 

IF!TIME .GT. 24.) TIME=!TIME - 24.)*2. + 24. 

IF!TIM£ .GT. 48.) TIME=!TIME - 48.)*3. + 48. 

T = TIME*3600. 

IF !TIME .GE. T23) GO TO 2 
DELT!K) = !ZU*T)**.5 

R 2= 1.45*! ! !OELTA*GRAV*VOL**2*T**1.5)/AA ) **.166667 ) 

HT = V0L/R2**2 
R2 = R2/100. 

RADIUS!K) = R2 

IF! RADIUS!K) .GE. RENO ) GO TO 9 
GO TO 21 

C SPREAD BY SURFACE TENS IGN-VI SCOUS 
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2 65= ZU*RHCIL*^2 

DtLT(K) = { ZU’f'T) *=i=.5 

R3 = ((SIGMA^^PZ =>= T**3 )/ BB )**.25 ^ 2.20 

HT = V0L/R3=«'=!'2 
R2 = R3/10Q. 

RACIUS(K) = R3 

IF( RADIUS(K) ,G6. REND ) GO TO 9 
GO TO 21 

9 RflCIUS(K) = REND 

HT = V0L/(REND*100.)*^‘2 
C 

C LOSS CUE TO EVAPORATION 

21 FF = C1=!'TIME=«'60./ ( 1. + C2’!‘TIME=!'60. ) 

VOL = { VVL - VVL*FF*0.01) 

VXL(K) = VOL/1000. 

C 

C PRINTS THE RADIUS OF OIL SPREAD IN NETERS, THE TIME OF 
C SPREAD IN HOURS, THE THICKNESS OF OIL AS A CUNCTICN OF 
C TINE IN CM, THE THICKNESS OF OIL LAYER IN CM, AND THE RE- 
C MAINING VOLUME OF OIL AFTER EVAPORATICN IN LITERS. 

HTT(K) = HT 

WRITE (6,103) RAOIUS(K) ,TIME ,HTT(K) ,0£LT(K) ,VXL(K) ,K 
X(K) = ALOGIO(TIME) 

Y(K) = AL0G10(RADIUS(K) ) 

ALPHA = 0.0 
TINES(K+1) = TIME 

C THE ESTUARINE NCNTIOAL DRIFT WAS SIMULATED AS A UNIFORM 
C CURRENT WITH CONSTANT SPEED APPROX. 65 M/HOURS IN THE 
C SOUTH BAY AND 200 M/HOURS IN THE NORTHERN BAY AREA. 

CRAN = 62.500«(TIMES(K+1) - TIMES(K)) 

IF( XMV(I) .GE. 16125. .AND. YMV(I) .GE. 51600. ) 
ICRAN = 208.3*(TIM£S(K+1) - TIMES(K)) 

1#(TIMES (K+1) - TIMES(K)) 

BETAl = 0.01745*( YFL^IOO. + 90. ) 

IF( XMV(K) .GT. 35475. .AND. YMV ( K ) .GT . 806 3 . 0 ) 

13ETA1 = 0.01745*( YFL^IOO. + 65. ) 

IF( XMV(K) .GT. 21736. .AND. YMV (K ) .GT. 17737. ) 

2BETA1 = 0.01745=«'( YFLY60. + 65. ) 

IF( XMV(K) .GT.21 736. .AND. YMV ( K ) .GT . 3 37 00 . ) 

38ETA1 = 0.01745^( YFL=^100. + 45. ) 

IF( XMV(K),GT. 16125. .AND. YM V ( K ) .GT .5 1 600 . ) 

4BETA1 = 0.01745*{ YFL't'60. + 160. ) 

IF( XMV(K) .LE. 16125. .AND. YMV ( K ) . LE . 51 60C . 0 ) 

5BETA1 = 0.01745*( YFL^130. + 90. ) 

XCRAN = CRANYCOS (BETAl ) 

YCRAN = CRANES IN ( BETAl ) 

IX = lY 

CALL RANDUdX, lY ,YFL ) 

BETA =0 .01745*( YFL^GAMAl + GAMA2 ) 

IF( XMV(K) .LE. 16125. ) 6ETA= 0 . 0 1745=!= ( YFL ^1 80 . + 9 0 . ) 
IF( XMV(K) .GT. 16125. .AND. YMV(K) .GT . 36230. ) 
IBETA = 0.01745*( YFL*GAMA3 + GAMA4 ) 

IF( XMV(K) .GT. 16125. .AND. YMV(K) .GT. 51600. ) 
2BETA = 0.01745^( YFL=i‘GAMA5 + GAMA6 ) 

ROUTE = 0SLV(K)*(TIMES(K+1)-TIM£S(K) ) 

XMV(K + 1) = XMV(K) + ROUTE=^COS (BETA) + XCRAN 
YMV(K + 1) = YMV(K) + ROUTE=!= S I N ( B ET A ) + YCRAN 
IX = lY 

C GENERATE RADIUS OF OIL SPREAD. 

DC 30 L = 1,KX 

ALPHA = ALPHA + 360 . / FLO AT ( KX ) 

CC(K,L) = XMV(K+1) + RADIUS (K)=!=CGS( ALPHA ) 

DC(K,L) = YMV(K + 1) + RADI US ( K ) *S I N ( ALPHA ) 

30 CONTINUE 

IF( HTT(K) .LT. 0.010 ) GO TO 22 
20 CCNTINUE 

22 WRITE(6,112) 

WPITE(6,202) 

D0 50 K=1,NK 

WRITE(6 , 201) WINOV(K) , OSLV( K) , SC WV(K) , XMV( K+1) , 
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2Y^V( K+1 ) tK 
50 CONTINUE 

V^PITE(6,112) 

C PICT THE SLICK MOVEMENTS INDUCED BY THE WIND 4NC ITS 
C SPREAC. 

NN = NK^KX 
L = 0 

DC 999 11=1, NK 
OG 999 JJ=1,KX 
L = L + L 

C LET CONVERT 50000 METERS EQUAL TO 10 INCHES: 

CCA(L) = CC( II, JJ ) 

CCAA(L) = 0.00020*(CCA(L) ) 

IF( CCAA(L) .LE. 0.0 ) CCAA(L) = 0.0 

IF( CCAA(L) .GE. 9.0 ) CCAA(L) = 9.0 

CCA(L) = DD( I I, JJ ) 

DCAA(L) = 0.00020=!'DDA(L ) 

IF( DDAA(L) .LE. 0.0 ) OOAA(L) = 0.0 

IF( CCAA(L) .GE.14.0 ) DOAA(L) = lA.O 

999 CONTINUE 
C 

C PLOT 

C 

C PLCT IN A LIN-LIN SCALE THE THICKNESS CF OIL AS A FUNCTION 
C OF TIME 

CALL PL0TP(X,HTT,NK,0) 

WRITE( 6, 112) 

C PLCT IN A LIN-LIN SCALE VOLUME AS A FUNCTION CF TIME. 

CALL PLCTP( X,VXL,NK,0 ) 

V<RITE(6,1121 

C PLCT IN A LOG-LOG SCALE THE RADIUS OF SPILL AS A FUNCTION 
C OF TIME 

CALL FL0TP(X,Y,NK,0) 

C 

REAL XTITLE( 2 )/' LATI ' , 'TUDE '/ 

REAL YTITLE( 3) / ' LONG' , • ITUD' , ' £ ' /’ 

XI = 0.0 
X2 = 0.0 
NCX = -8 
SIZEX = 8. 

XMIN = 0.0 
CX = 5C00. 

Y1 = 0.0 
Y2 = 0.0 
NCY = 9 
SIZEY = 14. 

YMIN = 0.0 
OY = 5000. 

DRAW 



C 



NK = NK+1 
CALL PLCTS 

CALL AXIS( XI, X2, XTITLE, NCX, SIZEX, 0.0, XMIN 
CALL AXIS( Y1,Y2, YTITLE,NCY, SIZEY, 90. , YMIN, 
NZ = 0 

CO 911 1=1, NK 

HTT(I) = O.QOOZO^i XMV(I) ) 

C£LT(I)= 0.0002-YMV(I ) 

NZ = NZ+1 

IF( HTT(I) .LT. 0.0 .OR. OELT(I) .LT. 0.0 
IF( HTT ( I ) .GE.8.0 .OR. DELT ( I ) . GE . 14.0 ) GG 
911 CONTINUE 

7 CALL LIN£(HTT,D£LT,NZ, 1, 2 ) 

CALL LINE (CCAA,DDAA,NN, 1,-7) 

CALL PLCT (0.0, 16. 0,-3 ) 



t DX) 
CY ) 



) GO TO 7 
TO 7 



IF( YMV(l) .GE. 36757.5 ) GO TO 555 
IF( XMV(l) .GE. 48165. ) GO tq 556 
IF( YMV(l) .GE. 60850. ) GO TO 557 
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IF( XMVd ) .GT. 58000. .OR. YMV ( 1 ) .GT. 7C900. J STOP 

SOUTHERN BAY 

DC 912 I=ltKK 
XC(I ) = X0(I ) - 140. 

!F( XQ( I ) .L£. 0.0 ) XQ(I ) = 0.0 
XG(I) = XQ(I )4=0.0725 
YQ(I ) = YOa ) 

IF( YQ(I) .GE. 125. ) YQ(I) = 125. 

YQ(I) = YQ(I )=i'0.0725 

912 CONTINUE 
SUM = 0.0 
NZ = 0 

CO 401 I=1,MK 

HTT(I) = 0.00039^XMV(I ) - 14.0 
CELT(I) = 0.00039*YMV( I ) 

IF(HTT(I).LT.O.O. OR . DELT ( I ) . LT . 0 .0 ) GO TO 401 
IF (HTT( I ) .GT .9.0 .OR. OELTd ) .GT. 12.0) GO TO 401 
CZd) = ( YMV(I) - YMV(I-l) )^^2 
HZ(I) = ( XMVd) - XMV( I-l ) )*=!=2 
NZ = NZ+1 

FA = 2.*RA0IUS(I j^^'SQRTi DZd) + HZ(I) ) 

SUM = SUM + FA 

IF( CELT(I) .LT. 3.5 ) GO TO 666 
C CONVERT TO THESIS FORMAT 8.5X11. 

HTTd) = HTT( 1)^0.725 
CELTd) = OELTd )*0.725 

IF( HTT ( I ) .GE.3.0 .OR. OELTd). GE. 8.4 ) GO TO 3 
401 CONTINUE 

8 IF( NZ .LT. 3 ) GO TO 888 

CALL LIN£( HTT,0£LT,NZ, 1, 2 ) 

NM = NZ+KX 
DC 501 I=1,NM 

CCAA(I) = 0.00039=i'CCA( I ) - 14. 

IF( CCAA(I) .LE. 0.0 ) CCAAd) = 0.0 
IF( CCAA(I) .GT. 4.0 ) CCAAd) = 4.0 
OCAA( I) = 0.00039’!‘DOA( I ) 

IF( DOAA(I) .GE.12.0 ) ODAA(I) = 12.0 
C CONVERT TO THESIS FORMAT 8.5X11. 

CCAAd ) = CCAA ( I ) *0.725 
. CCAA( I) = CDAA( I )*0.725 
501 CONTINUE 

SLMl = 0.0001*SUM 
WRITE(6t304) SUMl 
WRITE(6,303) RAOIUS(NZ), NZ 
CALL LINE( XO, YQ,KK, 1,1) 

CALL SYMBOL( 2.3, 9. 4, . 14, » EAST SOUTHERN BAY', .0,17 ) 
CALL LINE(CCAA,OOAA,NM, 1,-7 ) 

CALL PLOT (0.0, 12. 0,-3 ) 

C 

838 DC 913 1=1, KX 

XS( I ) = XOd ) - 60. 

IF ( XS( I ) .LE. 0.0 ) XS( I ) = O.C 
IF ( XS( I ) .GE. 90. ) XSd ) = 90. 

XSd ) = XSd )*0.0725 
YSd ) = YOd ) - 50. 

IF( YSd ) .LE. 0.0 ) YSd ) = 0.0 
IF( YSd) .GE. 120. ) YSd) = 120. 

YSd) = YSd )*0.0725 

913 CONTINUE 
SUM = 0.0 

HTTd) = 0.00039*XMVd) - 6.0 
DELT(l) = 0.00039*YMV( 1) - 5.0 

IF( HTTd).Gc. 0.0 .AND. OELTd). GE. 0.0 ) NZ = 0 
CO 701 1=1, NK 

HTTd) = 0.00039*XMV( I ) - 6.0 
OELTd) = C.00039*YMVd ) - 5.0 

IF (HTT (I ).LT.O.O. OR . CELT (I ) . LT . 0 . 0 ) GO TC 701 
IF(HTT( I ).GT.9.0. OR. CELT ( I ) .GT . 12.0 ) GO TO 701 
CZd) = ( YMV(I) - YMVd-l) )**2 
HZd) = ( XMVd) - XMVd-l))**2 



56 



I 




oooo 



NZ = NZ+l 

F4 = 2.*RA0IUS(I )*SQRT( DZ ( ! ) + HZ(I) > 

SUM = SUM + FA 

IF( HTT(I ) .GT.5.25.AiN0.DELT(I ) .GT.9.0 ) DEL7(I) 

IF( HTT(I),GT,8.0.AN0.DELT(I ).GT.6.5 ) OELT(I)= 

C CONVERT TO THESIS FORMAT 8.5X11. 

FTT(I) = HTT(I)«0.725 
DELT( I ) = D£LT( I )*0.725 
IF{ HTT(I) .LT. 2.4 .ANO.DELT(I) .LT. 7.3 ) GO TO 3 
IF(HTT( I).GE,5.5 .AND. OELT ( I ) .GE . 5. 0 ) GO TO 3 
701 CONTINUE 

3 CALL LINF(HTT,DELTtNZt It 2 ) 

NM = NZ^KX 
CO 601 I=1,NM 

CCAA(I) = 0.0C03 9=!'CCA( I ) - 6. 

IF( CCAA(I) .L£. 0.0 ) CCAA(I) = 0.0 
IF( CCAA(I) .GE. 9.0 ) CCAA(I) = 9.0 
CCAA(I) = 0.00039*D0A( I ) - 5.0 
IF( CCAA(I) .GE. 12.0 ) OOAA( I ) = 12.0 
IF( CCAA( n.GT.5.25.AN0.D0AA(I ) .GT. 9.0) DDAA(I)= 9.C 
IF( CCAA( I ).GT.3.0.ANO.DOAA(I) .GT.6.5 ) OOAA(I)= 6.5 
C CONVERT TO THESIS FORMAT 8.5X11. 

CCAAd ) = CCAAd )*0.725 
ODAA( I ) = DOAAd ) *0.725 
601 CONTINUE 

SLM2 = 0.0001*SUM 
WRITE(6»304) SUM2 
WPITE(6t303) RADIUS(iMZ), NZ 

CALL SYMBOL( 2.3, 9. 2 , . 14, • WEST SOUTHERN BAY', .0,17 ) 
CALL LINE( XS, YS,KK, 1, 1 ) 

CALL LINE(CCAA,DOAA,NM, 1,-7 ) 

CALL PLOT (0.0, 12. 0,-3) 



NORTHERN BAY 

556 DC 916 1=1, KK 

XVd) = XO(I) - 140. 

IF( XV( I ) .LE. 9.0 ) XVd ) = 0.0 
XVd) = 0.0725*XVd) 

YVd ) = YOd ) - 230. 

IF( YV( I ) .LE. 0.0 ) YVd ) =0.0 
YVd) = 0.0725*YVd) 

916 CONTINUE 
SUM = 0.0 
NZ = 0 

DC 705 1=1, NK 

HTTd) = 0.00039*XMVd ) - 14.0 
OELT(I) = 0. 00039*YMVd ) -23,0 
IF (HTT( I ) .LT .0.0 . GR . OELT (I ) .L T . 0 .0 ) GO TO 705 
DZd) = ( YMVd) - YMVd-1) )**2 
HZd) = ( XMVd) - XMVd-l))**2 
NZ = NZ+1 

FA = 2.TRADIUSd )*SQRT( DZd) + HZd) ) 

SUM = SUM + FA 
C CONVERT TO THESIS FORMAT 8.5X11. 

HTTd) = HTTd)*0.725 
DELT( I ) = OELTCI )*0.725 

IF( HTTd). GT.6.5 .OR. OELT d ). GT . 5.0 I GO TO 6 
705 CONTINUE 

6 IF( NZ .LT. 3 ) GO TO 666 

CALL LIN£( HTT,D£LT,NZ, 1, 2 ) 

NM = NZ*KX 
DO 602 1=1, NM 

CCAAd) = 0.00039*CCA(I ) - 14. 

IF( CCAAd) .LE. 0.0 ) CCAAd) = 0.0 

IF( CCAAd) .GE. 9.0 ) CCAAd) = 9.0 

DCAA(I) = 0.0C039*D0A d )- 23. 

IF( DCAA(I) .LE. 0.0 ) DDAAd) = 0.0 

IF( DDAA(I) .GS. 7.0 ) DDAA(I) = 7.0 

C CONVERT TO THESIS FORMAT 3.5X11. 

CCAAd) = CCAAd )*0.725 
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CC^A(!) = DDAA( I )=!'0.725 
602 CONTINUE 

SUMS = 0.0001*SUM 
WRITE(6,30A) SUMS 
WRITE(6,303) RADIUS(NZ), NZ 

CALL SYMBOL( 2 . 1 , S .43 , . 14 , • EAS T NORTHERN 8AV',.0, 
CALL LINE< XV,YV.KK, 1, 1 ) 

CALL LINE(CCAA,DDAA,NM, 1,-7 



CALL PLQT( 0.0, 10.0, -3 ) 



) 



XU( I ) 
XU( I ) 



0.0 

90. 



YU(I ) =0.0 



SS7 DO 91S 1=1, KK 

XU( I ) = XOd ) - SO. 

IF ( XU( I ) .LE. 0.0 ) 

IF( XU( I ) .GE. 90. ) 

XU(I) = 0.0725-XU(I) 

YU(I) = Y0(!) - 230. 

IF( YU( I ) .LE. 0.0 ) 

YUU) = 0.0725-YU(I) 

91S CONTINUE 
SUM = 0.0 

HTT(l) = 0.00039=f=XMV( 1) - S.O 
D£LT(L) = 0.00039*YMV( 1) - 23.0 
IF( HTT(1).GE. 0.0 .AND. DELT(l) 

DO 704 1=1, NK 

FTT(I) = 0.00039=i'XMV( I ) - S.O 
\ DELT(I) = 0.00039*YMV( I ) -25.0 
IF( HTT(I) .LE. 1.0 ) HTT(I)=1.0 
IF( OELT(I).G£. S.8 ) 0ELT(I)=S.8 
IF(HTT( I ) .LT.0.0 . OR. DELT ( I ) . LT . 0 .0 ) GO 

OR. OELT ( ! ) .GT . 6.0) GO 



GE. 0.0 ) NZ = 0 



6.3, .14, 'CENTRAL 

- YMV(I-l) )’i'’*'2 

- XMV( I-l) 



)*SQRT( DZd ) + hZ( I) ) 



IF(HTT( I ) .GT.9.0 . 

CALL SYMBOL( 2.2, 

DZd) = ( YMV(I) 

HZd) = ( XMV(I) 

NZ = NZ+1 
FA = 2. ’^RADIUSC I 
SUM = SUM + FA 
C CONVERT TO THESIS FORMAT 8.SX11. 
HTTd ) = HTT(I )*0.72S 
CELT( I ) = DELT( I )*0.72S 
IF( HTTd) .LE. 0.72. OR. OELT(I) 
704 CONTINUE 

S CALL LINE( H7T,DELT,NZ, 1, 2 ) 

NM = NZ*KX 
DO 604 I=1,NM 

CCAA(I) = 0.00039=«=CCA( I ) - S. 

IF( CCAA(I) .LE. 1.0 ) CCAAd) = 

IF( CCAA(I) .GE. 9.0 ) CCAAd) = 

OCAA( I) = 0.00039*DDA( I )- 23. 

IF( DDAA(I) .LE. 0.0 ) DDAA(I) = 

IF( DCAA(I) .GE. 5.3 ) DDAA(I) = 

C CONVERT TO THESIS FORMAT 8.5X11. 
CCAAd ) = CCAAd )=^0.72S 
CDAA( I) = DDAAd )’«'0.725 
604 CONTINUE 

SUM4 = 0.0001*SUM 
WRITE! 6,304) SUM4 
WRITE(6,503) RADIUS(NZ), 

CALL LINE! XU, YU,KK, 1, 



BAY' 



■0 704 
TO 704 
0,11 ) 



.GE. 3.75) GO TO 



l.C 

9.0 



0.0 

5.3 



NZ 
1 ) 



CALL SYMBOL! 2 . 1 , 5 .43 , . 14 , • WEST NORTHERN BAY'r.O, 
CALL LINE(CCaA,DDAA,NM, 1,-7 ) 

CALL PLOT! 0.0, 10.0, -3 ) 

CENTRAL BAY. 



555 DC 914 1=1, KK 

XT!I ) = XO! I ) - 40. 
IF! XT! I ) .LE. 0.0 ) 



XT!! ) = 0.0 



IF! 


XT! I ) 


.GE. 


90 . ) 


XT! I ) = 90. 


XT!I) = 0 


.0725*XT!I ) 




YT!I) = YO!I ) 


- 145. 




IF! 


YT! I ) 


.LE. 


0.0 ) 


YT! I ) =0.0 


IF! 


YT! I ) 


• GE • 


90. ) 


YT!I ) = 90. 



17 ) 



5 



17 ) 
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YTU» = 0.0725*Y7(I) 

914 CCNTINUE 
SUM = 0.0 

HTT(l) = 0.00039=!'XMV( 1) - 4.0 
OeLT(l) = 0.00039*YMV( 1) - 14.5 

IF( HTT(1).GE. O.O .ANO. DELKD.GE. 0.0 ) NZ = 0 
00 703 1=1, NK 

hTT(I) = 0.00039*XMV( I ) - 4.0 
OELT(I) = 0.00039*YMV( I ) -14.5 

IF(HTT( I ).LT.0.0. OR. DELT ( I ) . LT . 0 . 0 f GO TO 703 
IF(HTT( n.GT.9.0. OR. DELT ( I ) .GT . 9.0) GO TO 703 
CZ(I) = ( YMV(I) - YMV(I-l) 

HZU) = ( XMV(I) - XMV(I-1)3«*2 
NZ = NZ+1 

FA = 2.*RA0IUS(I )*SQRT( OZ(I) + HZ(I) ) 

SUM = SUM + FA 

IF(HTT( I ) .L7.2.125.AND.0ELT ( I ) . G7 . 2. 87 5 ) H7T( I )=2.125 
C C0NVER7 70 THESIS FORMAT 8.5X11. 

HTT(I) = hT7(I)*0.725 
OELT( I ) = OELTd )=f'0.725 

IF( H7T (I ) .GE.5.6 .OR. OELTd). GE. 6.5 ) GO TO 4 
703 CONTINUE 

4 CALL LINE{HTT,05LT,MZ, 1, 2 ) 

NM = NZ*KX 
DO 603 I=1,NM 

CCAA(I) = 0.00039*CCA( I ) - 4. 

IF( CCAA(I) .LE. 0.0 ) CCAA(I) = 0.0 
IF( CCAA(I) .GT. 7.3 ) CCAA(I) = 7.5 
CCAA(I) = 0.0C039*DDA( I )- 14.5 
IF( CCAA(I) .LE. 0.0 ) DOAA(I) = 0.0 

I F(CCAA( I ) .LT.2. 12 5.AN0.0DAA( I ) .GT.2.875 ) CCAA( I ) = 2.12 
IF( CCAA(I) .GE. 9.0 ) ODAA(I) = 9.0 
C CONVERT TO THESIS FORMAT 8.5X11. 

CCAA( I ) = CCAAd )=«'0.725 
DCAA( I) = DDAA( I ) *0.725 
603 CONTINUE 

SUM3 = 0.0001*SUM 
WRIT5( 6,304) SUM3 
WRITE(6,303) RADIUS(NZ), NZ 
CALL LINE( XT,YT,KK, 1, 1 ) 

CALL LIN£(CCAA,DOAA,NM, 1,-7 ) 

CALL PL0T(0. 0,12. 0,-31 



666 00 10 1=1, KK 

XPd ) = .0306*XO( I ) 

YPd ) = .0306*Y0 ( I ) 

10 CONTINUE 

DC 667 J=1,N 
AX(J) = AX(J)*0.78 
AY( J) =AY( J )*0.78 
667 CONTINUE 

00 668 K=1,M 
XR(K) = XR(K)*0.78 
YR(K) = YR(K)*0.78 
668 CONTINUE 

CALL LINE( XP, YP, KK ,-1,1) 

CALL SYMB0LC1.44, 10.0,.14, 

I'ANALYSIS OF OIL SPILL MOVEMENT ', .0, 30 ) 

CALL SYM60L(1.75, 9. 7, .14, 'IN SAN FRANCISCO BAY AREA', 
2.0,25 1 

CALL SYMB0L(0.77, 9. 4, .14, 

I'AS PREDICTED BY ESTUARINE NON TIDAL OR I FT ' ,0.0, 41 ) 

CALL SYMBOL(5.8, 7 . 95 , 0 . 07 , ' SUI SUN BAY', 0.0, 10) 

CALL SYMBGL(0.20, 7. 74 , . 07 , ' 38N ' , 0 .0 , 3 ) 

CALL SYMBOL( 2.30, 7 .78, .07, ' SAN OABLO BAY', 0.0, 13) 

CALL SYMBQL(5.06, 7. 54, .07,'* BEN IC I A ' , 0 .0 , 9 ) 

CALL SYMB0L(3.70,5.60,.07, '* A L B ANY ' ,0 . 0 , 8 ) 

CALL SYMBOL (4. 40, 4. 10 ,.14, 'ALAMEDA' ,0.0,7) 

CALL SYMBOL(2.17,3.73,0.07, 'HUNTER PC I NT ' , 0 . 0 , 1 2 ) 

CALL SYMBOL{3.5,2.72,.07, 'SAN FRANCISCO BAY', 0.0, 17) 
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CALL SYMBOLO. 50, 1.75 ,.07, 'FOSTER CITY «',.0,13) 
CALL SYMBOL (4,55 ,0.4,0.14, ' PALO ALTO', 0.0, 9) 

CALL L!N5(AX,AY, N, 1,1) 

CALL LINE(XR,YR, M, 1,1) 

CALL PLOKO. 0,14. 0,-3) 

CALL PLGTE 



STCP 
END 

//GO.SYSIN DO ’i' 



/ 
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COMPUTER SIMULATION PROGRAM OF 
THE SOCIAL COST OF THE LOSSES IN THE RECREATION 
ACTIVITIES AND THE COMMERCIAL FISHING ACTIVITIES 



//HAN1874 JOB (2696, 1446, RL44) , 'iMUDJIARCJO SMC 1874* 

// EXEC FORTCLGP, REGION. G0=180K 
//FGRT.SYSIN DO ’i' 

DIMENSION XD(50), YD(50), P(50), 7X ( 10 ) , TY ( 10 ) 
DIMENSION XH(50>, YH(50), XC(50», YC ( 50 I , YCC(50) 
DIMENSION XE(IO), Y5(10), PYY(IO), V(10), SC(IO) 

C INITIAL INPUT DATA 

REA0(5,89) ALFA, BEFA,PX,PY, AK, BK 
K = 4 

C READ VOLUME SPILLED IN LITERS. 

READ(5,83) 1V(J), J=1,K) 

WRITE(6,85 ) 

WRITE(6,82) (V(J), J=1,K) 

C 

79 FORMAT!/////) 

80 FORMAT!/, 10X,'N0: (• ,13, •)• ,2X, 'SOCIAL COST=', F12.2, 
12X,» VOL. SPILLED LTRS=',F12.1 ) 

81 FORMAT! / ,58X, • AH=' ,F9 .2, 3X, ' BH = ' , F9.2) 

32 FORMAT! /,10X,5F12.1 ) 

83 FORMAT! 5F12.1 ) 

84 FORMAT!/, lOX, 'PYY! ', 13,' ) SO. PRICE LEVEL OF Y:',F8.2) 

85 FORMAT! /,4X, ' INPUT DATA: ' ) 

36 FORMAT! » 1' ,4X, 'COORD. P~S OF PRGD.FNC, MAX. PROFIT', 

1' LINE, optimal PROFIT LINE AND WELFARE FUNCTION' ) 

87 FORMAT!//, 4X,' EQUILIBRIUM POINT: ' ) 

88 FORMAT !5X, ' AL PHA = ' , F4. 2 , 3X , ' BET A= ' , F4. 2 , 3 X , ' PX= ' , 
1F8.2,3X, 'PY=* ,F3.2,3X,' AK=' ,F9.2,3X, '8K=' ,F9.2 ) 

89 FORMAT! 2F4. 2, 2F 8. 2, 2F9.2 ) 

90 FORMAT!/, 10X,F9.1,3X,F9.1,5X,F9*1,5X,F9.1,3X,F12.1 ) 

91 FORMAT!/, 5X, '$ PROFIT OF X £ Y= ',F12.2,3X, 

I'WELFARE INITIAL STATUS :', F8. 1 ) 

92 FORMAT! / , lOX , F 3 . 2 , 3 X , 4F9 . 1 ) 

93 FORMAT !/, lOX, 'X UN ITS : ' , F9. 1 , 3 X , ' Y UNITS :', F9. 1 , 3X , 
l'J=',I3 ) 

94 FORMAT!'!',’ COORDINATE POINTS OF SUPPLY AND DEMAND', 
2' CURVE' ) 

N = 25 

M = 20 

CW = 2.50 

XH!11 = 0.0 

ZETA = ALFA + BEFA 

DO 1004 J=1,K 

C PRCDUCTION FUNCTION AS A FUNCTION OF VOLUME SPILLED. 

AKl =!V! J) /10.=<=*6 ) **2 

AH = AK*!0.9*£XP ! -AKl )+ 0.1) 

BKl = !5 .«V( J )/10 .*^6)=i'*2 

BH = BK*!0,75*EXP(-BK1) + 0.25) 

WRITE!6,79) 

WRITE(6,88 ) ALFA, BEFA, PX,PY,AK,BK 
WRITE(6,81) AH, BH 
XUP = PX*AH>i'=S‘2 
YUP = PY*BH**2 

QP = SQRT !PX**2-AH-'!==i‘2 + PY=<'=<'2’i=BH*=i'2) 

RL = XUP/QP 
RM = YUP/QP 
C EQUILIBRIUM POINT. 

XE(J) = SORT (ALFA^i'AH^ + Z/ ZETA ) 

YE!J) = SQRT!BEFA*8H=^*2/ ZETA ) 

C MAXIMUM PROFIT 

PROFT= PX=i'!XUP/OP ) + PY*(YUP/gP) 
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C WELFARE STATUS VALUE 

WWX = ALFA*PROFT/ ( PX=!=ZETA) 

WWY = BEFA*PROFT/ ( PY’i'ZETA) 

WS = CW^XE( J )*=!=ALFA YE(J)**B£FA 
WRITE(6,87) 

WRITE(6,93) RL,RM, J 
WRITE(6,93) WWX,WWY, J 
WRITE(6,93) XE(J), YE(J), J 
WRITE(6,91) PPOFT,WS 

C M3VING OF PRICE TO EQUILIBRIUM STATE. 

6UL = BH**2 - YE ( J 

IF( BUL.LT. 0.0) BUL=0.0 

PYY( J) = YE( J)*PX*AH/( BH=i'SQRT(BUL ) ) 

WRITE(6,84) J,PYY(J) 

WRITc(6,86) 

DO 1000 1=1, N 

C PPCDUCTICN POSSIBILITIES CURVE 

QRS = BH**2 - ( BH^=^2/AH*’!'2 )*XH( I )=!'=«'2 
IF( QRS .LT. 0.0) QRS = 0.0 
YH4I) = SQRT(QRS) 

C BUDGET CCMSTRAINT OR PROFIT FUNCTION 
YC(I) = (PROFT - PX*XH(i))/PY 
IF{ YC(I) .LT. 0.0 ) YC(I)=0.0 

C NEW OPTIMAL PROFIT LINE 

PROFC = PX=!'XE ( J ) + PYY(J)^YE(J) 

YCCd )=(PROFC - PX=i'XH( I ) )/PYY( J ) 

IF( YCCd ) .LT. 0.0) YCCd)=0.0 

C WELFARE CURVE 

Xhd + 1) = XHd) + 50. 

XOd ) = XHd ) 

IF( XD(I) .EQ. 0.0 ) XD(I) = 10.0 
YCd) = (WS/(CW^=XD( I )=^=*ALFA ) )**< 1./3EFA ) 

WRITE (6, 90 ) XH( I ) , YH( I ) ,YC ( I ) , YCC (I ) ,YD( I ) 

1000 CONTINUE 
WRITE(6,94 ) 

C 

C CONVERT TO THESIS FORMAT 
DO 1001 I=1,N 

Xh(I) = 0.00005*XH(I ) + 0.5 
IF( XH( I ) .GT .6.5 ) XH(I)=6.5 
YFd) = 0.00005=i‘YHd ) + 0.5 
IF( YH( I ) .GT.9.0) YH(I)=9.0 
YCd) = 0.0Q005=!'YC (I ) + 0.5 
IF( YC(n. GT.9.0) YC(I)=9.0 
XCd) = 0.00005*XD(I) + 0.5 
IF( X0( I ) .GT.6.5 ) XOd )=6.5 
YDd ) = 0.005=!'YDn ) + 0.5 
IF( YD( I ) .GT.9.0) YD(I)=9.0 
YCCd)= 0.00005=i'YCCd )+ 0.5 
IF( YCC ( I ) .GT.9.0 ) YCCd) = 9.0 

1001 CONTINUE 
TXd) = 0.5 
TY(1) = 0.5 
TX(2) = 6.5 
TY(2) = 0.5 
TX(3) = 6.5 
TY(3) = 9.0 
TX(4) = 0.5 
TY(4) = 9.0 
TX(5) = 0.5 
TY(5) = 0.5 
CALL PLOTS 

CALL SYMBOLC 2 . 5 , 0. 0 , . 14, ' X UNIT PRODUCED ',0.0,17 ) 
CALL SYMBOL! 0 .0 , 2 .0 , . 14, • Y UNIT PRODUCED ’,90.0,17 ) 
CALL LINE! TX,TY,5,1,1 ) 

CALL SYMBOL!4.0,8.0,.07, PROD. PCSS CURVE ’,0.,19) 
CALL SYMB0L!4.0,7.3,.07,’- MAX PROFIT LINE ',0.0,18) 
CALL SYMBOL! 4.0, 7.6, .07, WELFARE CURVE ',0.0,16) 

CALL LINE! XH,YH, N,l,2 ) 

CALL LINE! XH,YC, N,l,-7 ) 

CALL LINE! XH,YCC,N,1,6 ) 

CALL LINE! XD,YD, N,l,5 ) 
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CALL PLCT{ 0.0, 12. 0,-3 ) 



P(l) = 5. 

DC 1002 I=1,M 
SUPPLY CURVE 

XH(I) = AH*=!=2*P( I ) /PROFC 
YH(I) = 6H*=i'2*P{ I ) /PRCFC 
DEMAND CURVE 

XC(I) = ALFA«PR0FC/ZETA*( l./Pd ) ) 

YC(I) = B£FA=!'PR0FC/ZETA=i=( l./P( I ) ) 

WRITE (6 ,92) P( I ) , XH( I ) ,YH( I ) ,XC (I ) ,YC( I ) 

P(I+1) = P(I) + 5. 

1002 CCNTINUE 

CONVERT TO THESIS FORMAT 
DC 1003 I=1,M 
PCI ) = 0.025*P( I ) + 0.5 
IF( P( I ).GT.9.0) P(I )=9.0 
Xhd) = 0.00005*Xh(I ) + 0.5 
IFCXHCI ).GT.6.5) XH(I) = 6.5 
XC( I ) = 0.00005*XC( I ) +0.5 
IFCXCd ).GT.6.5) XCd ) = 6.5 
YCd) = 0.00005*YC (I ) + 0.5 
IFCYCd ).GT.6.5) YCd ) = 6.5 
Yhd) = 0.00C05*YHd ) + 0.5 
IF(YH( I ).GT.6.5> YH(I ) = 6.5 

1003 CONTINUE 

CALL SYMBOL! 2 .5 ,0.0 , . 14, * X UNIT PRODUCED ',0.0,17 
CALL SYMBOL! 0 . 0 , 2 .0 , . 14, • PRICE PER UNIT PRODUCT ' 
190.0,24 ) 

CALL SYMBOL!4.0,3.2,.07, ' 

CALL SYMB0L!4. 0,8.0,. 07, ' 

CALL LINE! XC , F , M,l,l ) 

CALL LINE! XH,P,M,1,2 ) 

CALL LINE! TX,TY,5,1,1 ) 

CALL PLOT! 0.0, 12. 0,-3 ) 



SUPPLY CURVE • ,0.0, 16 ) 
DEMAND CURVE • ,0.0, 16 ) 



CALL SYMBOL! 2 .5 ,0 .0, . 14, • 
CALL SYMBOL! 0 .0 ,2 .0 , . 14, ' 
190.0,24 ) 

CALL SYMBJL!4.0,8.2,.07, ' 
CALL SYMB0L!4.0,8.0,.07, • 
CALL LINE! YC, P, M,l, 1 ) 
CALL LINE! YH,P, M,l,2 ) 
CALL LINE! TX,TY,5,1,1 ) 
CALL PLOT! 0.0, 12. 0,-3 ) 
CALL PLCTE 
1004 CCNTINUE 



Y UNIT PRODUCED • ,0.0, 17 
PRICE PER UNIT PRODUCT ' • 



+ SUPPLY CURVE 
- DEMAND CURVE 



' ,0.0,16) 
• ,0.0, 16 ) 



SOCIAL COST OF RECREATION ACTIVITIES C COMMERCIAL FISHING 
ACTIVITIES !AND RELATED INDUSTRIES) 

WRITE!6,79) 

DO 1005 J=1,K 

SC!J) = PX#!XE!1) - X£!J)) + P YY ! 1 ) ^ ! YE ! 1 ) - Y£!J)) 
WRIT£!6,80) J, SC!J), V!J) 

1005 CCNTINUE 
STOP 
END 

//GO.SYSIN DC * 
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